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AN  ASSESSMENT  OF  THE  USE  OF  STRAIN  AND  INERTIAL  SEISMOGRAPHS 
TO  ENHANCE  SEISMIC  SIGNAL  TO  NOISE  RATIOS 


1 .  FOREWORD 


In  order  to  verify  compliance  with  the  terms  of  the  proposed  comprehensive 
nuclear  test  ban  treaty,  seismic  monitoring  sites  may  be  located  at  regional 
distances  from  potential  test  sites.  It  can  be  reasonably  assumed  that  some 
type  of  signal-to-noise  ratio  (SNR)  enhancement  technique  wilj,  have  to  be 
applied  to  the  raw  data  collected  at  these  sites  in  order  to  detect  and  iden¬ 
tify  the  source  of  the  relatively  weak  seismic  signals  generated  by  low 
magnitude  earthquakes  and  underground  nuclear  explosions.  In  the  past, 
seismic  arrays  coupled  with  simple  beamforming  data  processing  methods  have 
been  used  extensively  to  enhance  short  period  teleseisraic  SNRs.  Utilization 
of  this  technique  is  based  upon  the  observation  that  in  the  frequency  band¬ 
width  commonly  occupied  by  teleseismic  signals,  the  signal  field  generally 
appears  to  be  much  better  organized  than  the  ambient  earth  noise.  However, 
recent  studies  by  Mrazek  et  al  (1979)  indicate  that  the  lateral  decorrelation 
of  regional  signal  fields  is  comparable  to  that  of  the  earth  noise  field  in 
commonly  shared  bandwidths.  Under  these  conditions,  the  use  of  horizontally 
distributed  arrays,  coupled  with  simple  beamforming  techniques  is  not  likely 
to  yield  significant  improvements  in  regional  SNRs.  Thus,  there  is  a  need 
to  reassess  the  use  of  existing  SNR  enhancement  methods  as  they  pertain  to 
the  detection  and  identification  of  weak  regional  seismic  signals.r^In  this 
regard  it  is  worthwhile  to  reconsider  the  use  of  linear  combinaticStvof  the 
outputs  of  strain  and  inertial  seismographs  to  increase  regional  ^NRs.  Like 
seismic  arrays  this  particular  method  utilizes  differences  in  apparent  hori¬ 
zontal  wavenumber  to  separate  signals  from  noise.  However,  in  contrast  to 
arrays,  it  is  insensitive  to  the  effects  of  lateral  signal  decorrelation 
because  all  necessary  observations  are  made  at  essentially  the  same  point. 
Colocation  of  the  sensors  may  also  substantially  reduce  the  potential  politi¬ 
cal,  logistical  and  security  problems  associated  with  the  operation  of  a 
nuclear  test  monitoring  station  in  a  foreign  country.  Thus,  the  strain- 
inertial  (S-I)  SNR  enhancement  method  seems  particularly  well  suited  for  use 
at  the  treaty  verification  sites  that  may  be  located  within  the  Soviet  Union. 
Unfortunately,  previous  attempts  to  utilize  S-I  methods  in  the  past  have  met 
with  only  limited  success  and  investigations  in  this  area  were  terminated  in 
the  early  1970s  (cf.  Shopland  1966,  1968,  Shopland  and  Kirklin  1969,  1970). 
This  relatively  poor  performance  has  been  attributed  primarily  to  the 
existence  of  high  system  noise  levels  at  the  output  of  the  strain  transducer. 
However,  in  retrospect,  the  failure  to  recognize  and  adequately  account  for 
the  effects  of  local  environmental  factors  may  have  played  an  equally  impor¬ 
tant  role  in  limiting  the  SNR  enhancement  capabilities  of  previously  devel¬ 
oped  S-I  systems.  Recent  progress  in  the  development  of  high  resolution 
seismic  sensors  suggests  that  it  may  now  be  possible  to  substantially  reduce 
the  strain  system  noise  problems  that  plagued  earlier  attempts  to  reduce  the 
S-I  method  to  practice.  Furthermore,  the  greater  appreciation  of  the  local 
environmental  noise  problem  that  has  developed  during  the  past  decade  leads 
to  the  belief  that  practical  methods  for  minimizing  their  Impact  on  S-I  SNR 
enhancement  capabilities  can  be  developed. 
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For  these  reasons  AFOSR  und^t  the  sponsorship  of  DARPA  has  authorized  Teledyne 
Geotech  to  assess  the  merit  of  developing  an  advanced  S-I  seismograph  system 
to  enhance  seismic  SNRs.  The  principal  results  of  this  assessment  are 
described  in  the  following  pages. 


2.  POTENTIAL  P  WAVE  SNR  ENHANCEMENT 


2.1  BASIC  CONCEPTS 

Interest  in  the  strain-inertial  SNR  enhancement  method  was  triggered  by  Romney 
(1964)  who  presented  theoretical  arguments  to  demonstrate  that  simple  subtrac¬ 
tion  of  an  appropriately  filtered  and  normalized  strain  seismogram  from  a 
simultaneously  recorded  inertial  seismogram  would  result  in  substantial 
improvements  in  P  wave  SNRs,  provided  that  the  earth  noise  consisted  princi¬ 
pally  of  single  mode  Rayleigh  waves.  The  essential  elements  of  his  arguments 
are  illustrated  by  the  two  sets  of  curves  shown  in  figures  1-a  and  1-b.  For 
the  purposes  of  this  illustration  we  have  assumed  that  vertical  strain  and 
displacement  observations  are  made  at  the  surface  of  a  homogeneous,  isotropic 
perfectly  elastic  half  space.  We  have  further  assumed  that  the  earth  noise 
consists  of  Rayleigh  waves  and  that  we  are  interested  in  enhancing  the  SNR 
for  a  P  wave  arriving  at  the  observation  point  with  a  20°  angle  of  incidence. 
The  two  curves  shown  in  figure  1-a  illustrate  the  relative  differences  between 
the  strain-displacement  response  ratios  for  the  Rayleigh  wave  noise  and  the  P 
wave  signals.  These  ratios  measure  the  magnitude  of  the  vertical  strain  per 
unit  of  vertical  displacement  and  thus  reflect  the  relative  strain  sensitivity 
of  the  earth  to  different  types  of  seismic  waves.  For  the  purposes  of  illus¬ 
tration  they  have  been  normalized  so  that  the  Rayleigh  wave  strain-displacement 
response  ratio  is  1  at  a  frequency  of  1  Hz.  Observe  that  the  relative  verti¬ 
cal  strains  for  the  Rayleigh  wave  are  almost  a  factor  of  10  larger  than  those 
for  the  P  wave.  This  means,  in  effect,  that  the  output  of  a  vertical  strain 
sensor  will  be  relatively  rich  in  Rayleigh  wave  noise  when  compared  to  the 
output  of  the  vertical  displacement  sensor.  In  other  words,  the  strain  out¬ 
put  is  for  all  practical  purposes,  an  estimator  of  the  noise  observed  at  the 
output  of  the  vertical  displacement  sensor.  We  may  take  advantage  of  this 
property  by  filtering  the  strain  output  so  that  it  is  identical  to  the 
displacement  output  when  only  noise  is  present.  When  this  is  done  the  strain 
displacement  response  ratios  will  take  the  form  shown  in  figure  1-b.  From 
this  illustration  it  can  be  seen  that  if  the  filtered  strain  output  is 
subtracted  from  the  displacement  output  the  remainder  will  contain  only  the 
signal  plus  a  relatively  weak  signal  distortion  term.  Such  a  process  is  com¬ 
monly  referred  to  as  prediction  error  filtering  and  under  appropriate 
background  noise  conditions  this  technique  can  be  a  powerful  tool  for  seismic 
SNR  enhancement. 

In  the  following  sections  we  shall  quantitatively  assess  the  potential  SNR 
enhancement  capabilities  of  a  linear  combination  of  the  outputs  of  a  vertical 
strain  sensor  and  a  vertical  displacement  sensor  under  certain  idealized  noise 
conditions.  The  results  of  this  assessment  will  not  only  provide  insight  into 
the  problems  encountered  in  previous  attempts  to  reduce  the  S-I  method  to  prac¬ 
tice  but  will  also  provide  basic  constraints  on  the  design  and  utilization  of 
an  advanced  S-I  system. 
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2.2  CONCEPTUAL  NOISE  MODELS 

Realistic  assessment  of  the  SNR  enhancement  capabilities  of  S-I  systems 
requires  quantitative  information  regarding  the  frequency  wave  number  struc¬ 
ture  and  composition  of  the  earth  noise  in  the  regional  signal  bandwidth. 
Unfortunately,  a  review  of  the  existing  data  base  earth  noise  (see  section  6) 
has  revealed  that  this  type  of  information  is  generally  lacking,  particularly 
in  the  regional  signal  bandwidth.  Thus  we  must  resort  to  generalized  concep¬ 
tual  models  of  the  earth  noise  in  order  to  obtain  a  preliminary  assessment  of 
the  SNR  enhancement  capabilities  of  S-I  systems.  For  the  purposes  of  this 
investigation  we  shall  confine  our  attention  to  propagating  earth  noise  fields 
which  consist  of  both  scattered  Rayleigh  waves  and  P  waves.  We  shall  assume 
that  these  fields  are  laterally  homogeneous  and  time  stationary.  In  addition, 
the  effects  of  locally  generated  wind  noise  and  strain  system  noise  will  be 
evaluated  under  special  circumstances.  Specific  cases  examined  during  the 
course  of  this  investigation  are  listed  in  table  I.  It  will  be  noted  that 
SNR  enhancement  and  noise  rejection  capabilities  are  considered  separately. 
This  is  because  the  relationship  between  SNR  enhancement  and  noise  rejection 
is  generally  complicated  by  the  presence  of  the  signal  distortion  term 
referred  to  earlier. 

For  the  purposes  of  this  project  we  have  also  limited  ourselves  to  the  case 
where  the  observation  point  is  located  at  the  surface  of  a  homogeneous, 
isotropic,  elastic  half  space.  While  this  assumption  influences  the  details 
of  our  calculations  it  has  little  impact  upon  the  general  nature  of  our 
results.  Furthermore,  the  current  experimental  data  base  concerning  earth 
noise  does  not  justify  the  use  of  more  sophisticated  earth  models. 


2.3  PREDICTION  ERROR  FILTER  DESIGN 

For  the  remainder  of  this  discussion  let  us  focus  our  attention  upon  the 
potential  noise  reduction  and  SNR  enhancement  capabilities  of  a  linear  com¬ 
bination  of  the  outputs  of  a  vertical  inertial  and  vertical  strain  sensor. 

Let  ne(t)  and  n^(t)  be  the  outputs  of  the  strain  and  inertial  sensors, 
respectively,  when  only  noise  is  present.  We  wish  to  apply  a  filter  to  the 
output  of  the  strain  sensor  such  that  it  is  the  optimum  predictor  of  the 
noise  observed  at  the  output  of  the  inertial  sensor  in  the  sense  that  the 
mean  square  difference  between  the  observed  and  predicted  inertial  outputs  is 
minimized.  If  H(w)  is  the  complex  frequency  response  of  this  filter  then 

..(»)  =  (1) 

♦eeCw) 

where  4>ei(ji>)  is  the  crosspower  spectral  relating  the  strain  and  inertial  out¬ 
puts  and  <f>ee(w)  is  the  power  spectral  density  of  the  strain  output  (For 
readers  who  are  unfamiliar  with  two  channel  prediction  error  operators,  eq.  1 
is  derived  in  detail  in  section  1  of  the  Appendix). 


It  is  convenient  for  our  purposes  to  define  the  power  and  crosspower  spectral 
density  functions  in  terms  of  their  frequency  wave  number  spectra.  Thus, 


Table  1.  Strain-Inertial  test  cases 

1.  Noise  rejection  capabilities 

•  Rayleigh  wave  noise  contaminated  by  strain  "system"  noise 

•  Rayleigh  wave  noise  contaminated  by  "wind"  noise 

•  Rayleigh  wave  noise  contaminated  by  "mantle"  P  wave  noise 

2.  SNR  enhancement  capabilities 

•  Scattered  Rayleigh  wave  noise 

•  Scattered  P  wave  noise 


$ei(w,Ar) 


1  yNel  (co,k)  exp  (-jk*Ar)dk 
4  IT2  K 


(2) 


where  Nej  (w,li)  is  the  3  dimensional  crosspower  spectral  density  function,  k 
is  the  vector  horizontal  wave  number,  Sr  is  the  vector  distance  between  the 
strain  and  inertial  observation  points  and  the  operation 


(3) 


indicates  aj^  integration  over  the  entire  k  plane.  Similarly,  if  Nee  (u>,k) 
and  (u>,k)  are  the  3  dimensional  power  spectral  density  functions  charac¬ 
terizing  the  vertical  strain  and  inertial  noise  fields,  then 

"  r  _»  — *  —A  — k 

$ee  (w,Ar)  =  1  J  Nee  (rn,k)  exp(-jk.Sr)  dk  (4) 

47?  * 


and 


n 

4> it  (w.Ar) 


=  1  j  (w,k)  exp(-jk*4r)  dk 

5  J K 
4n^ 


(5) 


2.3.1  Rayleigh  Wave  Noise  Prediction  Filters 

Let  Mg  (w)  and  (m)  be  the  complex  frequency  responses  of  the  strain  and 
inertial  sensors.  If  the  propagating  earth  noise  field  consists  entirely  of 
single  mode  Rayleigh  waves  from  remote  sources  then  using  the  results  of  sec¬ 
tion  4.2  in  the  Appendix  it  can  be  shown  that 

4>ei  -  HI"  Mj  G*  Me  JIG  JG  Jr)  (6) 

3Z  r 


where  J0  ( tc r 6 r )  is  a  Bessel  function  of  order  0,  |l|2  is  relate^to  the 
distribution  and  spectral  intensity  of  the  sources  and  G*  and  -gj  are  the 
displacement  and  strain  medium  response  functions.  Similarly  if>ye  and  i/>§i  are 
defined  by 


C.  -  IM2M2la£l2  <'> 

+ii  -  IlI?Im1I"Ig!/  (8) 


Therefore  from  equations  1 

H(w) 


6,  and  7 


Me 

7 


3G 

3Z 


I  Me 


3G 

~2 

3Z 

(k«r) 


(9) 


It  is  important  to  note  that  in  this  case  when  strain  and  inertial  sensors 
are  colocated,  H(io)  is  independent  of  the  intensity  and  distribution  of 
sources  and  the  state  of  organization  of  the  corresponding  Rayleigh  wave  fields. 
This  means  in  effect  then  that  when  the  earth  noise  consists  solely  of 
unimodal  Rayleigh  waves  the  noise  prediction  filter  is  dependent  only  on  the 
strain  and  displacement  medium  response  functions  and  frequency  response 
functions  of  the  filters.  Since  the  parameters  which  determine  these  functions 
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are  time  invariant,  the  Rayleigh  noise  prediction  filter  also  shares  this 
property.  In  the  special  case  where  the  strain  and  inertial  sensors  are 
colocated  at  the  surface  of  a  homogeneous  and  isotropic  halfspace  the  medium 
response  functions  are  defined  by  equations  3.3.5  in  section  3  of  the 
Appendix  and  H(w)  reduces  to 


H(<o) 


M|  Me 

TCP 


[i-^rj2]1^  [2 


_ xr _ 

'jLy-iMWlrPPP 

£rj  \6  ) 


where  cr  is  the  Rayleigh  wave  speed.  If  we  assume  that  X  *  p  then  a  =  3^6~ 

and  cr  =  0.91946  and 

★  . 

H(w)  =  Mj  Me  I  2.337 ct  j  (11) 

I  Me  | 2  1  «  i 

The  term  in  brackets  is  plotted  in  figure  2  for  the  case  where  a  =  4.0  km/s. 


2.3.2  P  Wave  Noise  Prediction  Filters 


Let  us  now  assume  that  the  propagating  noise  field  consists  entirely  of  scat¬ 
tered  P  waves  with  apparent  horizontal  velocities  c  distributed  in  the  inter¬ 
nal  c,^,^  c  <_  a  where 


Cmin  >  “ 


(12) 


Using  the  3  dimensional  power  and  cross  power  spectral  density  functions  given 
in  section  4.1  of  the  Appendix,  in  equations  2,  4  and  5  it  is  found  that 

n  _i  b  2t\ 

!>ei  (w,<Sr)  =  MP*M-j  f  j  Ai  (0,u>)  3D  *  D  exp  (-jk.6)  kdedk  (13) 

oo  9Z 

b  2  it 


♦ 


ee  (MO  =  I 1  ?  ffAj  (0,w) 

n  n 


0  0 
b  ?■  TT 


exp  (-jk*6r)  kdOdk 


t>ii  (w,  5r)  =  I  Mi  |  ^  f  f  AA  (e,u>)  |  D  |  2  exp  (-jk.fir)  kdOdk 

4P  J  J 


(14) 

(13) 


o  o 


where  Ai  (0,oj)  is  the  power  spectral  density  of  the  displacement  potential 
for  the  incident  P  wave  field,  D  ig^the  vertical  displacement  reflectivity 
function  for  incident  P  waves  and  the  vertical  strain  reflectivity 

function.  In  the  special  case  where  the  strain  and  inertial  sensors  are 
colocated  at  the  surface  of  a  homogeneous,  and  isotropic  half  space  the 
results  of  sections  3  and  4  of  the  Appendix  may  be  used  to  show  that 
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FIGURE  2.  MODULUS  OF  THE  FREQUENCY  OF  A  RAYLEIGH  WAVE  NOISE  PREDICTION  FILTER. 
STRAIN  AND  INERTIAL  SENSOR  RESPONSES  ARE  ASSUMED  TO  BE  IDENTICAL. 
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¥ 


(16) 


where 


D  -  2J(| 


3D 

3Z 


4*  k2 

A  +  2p  a 


n  =  sin  6 


k  -  “ 
a  a 


(l-n2)1/2^(g)  -  2ri2J 

4,2  dV)  >/*((«)  2-n2)  >/2^(aj  2_^y 

"2  '‘-"2>I/2((t)2-"2)1/2 


4n2  (1-n2)  ^2 


(17) 

(18) 
(19) 


and  0  is  the  angle  of  incidence  for  the  P  wave  noise.  Given  these  definitions 
and  assuming  that  the  sensors  are  colocated  equations  13-15  may  be  written  as 

<f>ei  -- 8j  MP%  Aj0^)  fa\ 2  A  Iei  (nQ)  (20) 

4it2  X+2u 


where 


(21) 

(22) 

(23) 

(24) 

(25) 

(26) 
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Observe  that  in  this  case  the  noise  prediction  filter  is  independent  of  the 
azimuthal  distribution  of  the  P  wave  field  but  is  functionally  dependent  on 
no.  Now  ri0  is  related  to  c  m^n  through 

r  .  «  (29) 


Thus  the  range  of  apparent  horizontal  phase  velocities  for  the  scattered  P 
waves  contained  in  the  field  is  given  by 


Ot  <  Q  <  00 

~  (30) 

Ih(w)|  for  representative  values  of  n  are  plotted  in  figure  3  for  the  case 
where  Mg^M^,  and  a=4.0  km/sec.  Notice  that  the  gain  of  the  frequency 

response  functions  diminish  as  n  increases.  In  other  words,  in  a  P  wave 
noise  field,  the  shape  and  phase  of  the  noise  prediction  filter  response  are 
constants  but  the  gain  will  change  in  response  to  changes  in  the  range  of 
apparent  horizontal  velocities  contained  in  the  field. 


2.3.3  Prediction  Filters  for  Multicomponent  Noise  Fields 

Let  us  consider  the  case  where  the  outputs  of  the  strain  and  inertial  sensors 
contains  two  statistically  independent  components,  i.e.; 


ne(t)  =  nel (t)  +  ne2(t) 

(31) 

ni(t)  -  ni5 (t)  +  nj?(t) 

(32) 

Then 

<f>  gi(u)  “  <t>gi(w)  +  <f>el(u>) 

(33) 

and 

4>"e(w)  =  ^(w)  +  <t>ee(w) 

(34) 

1  2 

where  and  are  crosspower  spectral  densities  relating  the 

inertial  outputs  for  components  1  and  2  and  ^^(w)  +  'f'geCw)  are 
power  spectral  density  functions.  Therefore, 

strain  and 
the  strain 

4> 1 .  +  4>?  < 

H(w)  -  el  ei 

’ec  ce 

■  B]Hj(ui)  +  B^H.(ai) 

(35) 

where 

4,1 

f  ’ 

ee 

(36) 

4>? 

IL(«‘>)  =  ~7“ 

<f>cc 

(37) 
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FIGURE  3.  MODULUS  OF  THE  FREQUENCY  RESPONSE  OF  APWAVE  NOISE  PREDICTION 
FILTER.  STRAIN  AND  INERTIAL  SENSOR  RESPONSES  ARE  ASSUMED  TO  BE 
IDENTICAL 

o  mu 

-li¬ 


ra  80-6 


are  Che  optimum  frequency  response  functions  that  would  apply  for  fields  con¬ 
sisting  solely  of  component  1  or  component  2  and 


(38) 


(39) 


are  the  ratios  of  the  power  in  components  1  and  2  relative  to  the  total  noise 
power  observed  at  the  output  of  the  strain  sensor.  By  analogy,  then  if  the 
noise  field  consists  say  £  statistically  independent  components 


<u 

r  -  ee 

1  '  T~ 

ee 


V 


4>2 

T 


ee 


L 

H(w)  “  £  B  H  (u)  (40) 

«,  =  ! 


That  is,  the  optimum  noise  prediction  filter  for  a  multicomponent  field  is 
the  weighted  sum  of  the  optimum  prediction  filters  for  the  individual  com¬ 
ponents  where  the  weights  are  the  fractional  noise  power  in  each  component. 


Now  suppose  one  of  the  components  is  strain  system  noise,  then  for  the  case 
of  the  two  component  noise  field  ) 


H(w)  -  Bj  Hj  (w)  -  (1-BS)  Hj  (u) 


(41) 


since  the  strain  system  noise  has  no  counter  part  at  the  output  of  the  iner¬ 
tial  sensor.  Thus  the  gain  of  the  optimum  filter  will  decrease  as  the  per¬ 
centage  of  system  noise  in  the  total  strain  noise  increases. 


2.4  NOISE  REDI  CTION 


Let  be  the  power  spectral  density  of  the  output  of  the  noise  prediction 
filter.  It  follows  from  the  previous  discussion  that 

♦ii  -  Ih2I  *ee  (42) 


Le_  be  the  ratio  of  predicted  inertial  noise  power  to  the  actual  inertial 
noise  power.  Then 


Yei 


♦  ii 
<p  ii 


4> 


ii 


(43) 


The  terra  on  the  RHS  of  equation  43  is  recognized  as  the  square  of  the  simple 
coherence  relating  the  strain  and  Inertial  outputs.  Now  it  is  shown  in  sec¬ 
tion  1  of  the  Appendix  that  the  power  spectral  density  of  the  noise  remaining 
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after  subtracting  the  predicted  from  the  observed  inertial  output  is  given  by 

N“*ii(l-Yei>  <44> 

The  ratio  of  the  noise  power  remaining  after  prediction  error  filter  to  the 
observed  inertial  noise  power  is  therefore 

R  -  N  -  i-y2ei  (45) 

7T  ei 

R  measures  the  noise  reducing  capabilities  of  the  prediction  error  filtering 
process.  In  the  following  paragraphs  we  shall  investigate  the  behavior  of  R 
under  various  seismic  noise  conditions. 

2.4.1  Rayleigh  Waves  Contaminated  by  Strain  System  Noise 

It  is  shown  in  section  5.1  of  the  Appendix  that  if  the  earth  noise  consists 
entirely  of  Rayleigh  waves  of  the  same  mode  then 

R  *  B  (46) 

where  B  is  the  ratio  of  system  noise  power  to  total  noise  power  observed  at 
the  output  of  the  strain  sensor.  It  is  important  to  notice  that  in  this  case 
the  noise  reducing  capabilities  are  inherently  broad  band  and  limited  only  by 
the  presence  of  system  noise.  Clearly  if  strain  system  noise  levels  can  be 
maintained  at  relatively  low  levels  then  a  linear  combination  of  vertical 
strain  and  inertial  transducers  can  be  a  powerful  tool  for  reducing  Rayleigh 
wave  noise. 

2.4.2  Rayleigh  Waves  Contaminated  by  Wind  Noise 

Earth  movement  in  response  to  wind  generated  atmospheric  pressure  changes  is 
a  common  constituent  of  the  earth  noise  observed  at  many  surface  locations. 

It  is  worthwhile  therefore  to  determine  its  impact  on  the  noise  reducing 
capabilities  of  a  vertically  oriented  S-I  system.  It  is  shown  in  section  5.2 
of  the  Appendix  that  in  this  instance 


( 

\  1-T 

r  b(i-b) 

(-?) 

2 

B(l-B)  + 

V| 

(i-t)b 

where  B  is  the  ratio  of  the  wind  noise  power  to  the  total  noise  power  at  the 
output  of  the  strain  sensor  and  Q-1  and  T-*  are  the  moduli  of  the  strain- 
displacement  response  ratios  for  the  wind  noise  and  Rayleigh  wave  noise  com¬ 
ponents  respectively.  It  is  also  shown  in  section  5.2  that  if  the  observation 
point  is  at  the  surface  of  a  homogeneous  and  isotropic  half  space  then  £  is  of 
the  order  of  10”^  and  that  T 


R  =  B  ,  B  t  1  (48) 

Equation  47  has  been  evaluated  explicitly  for  the  case  where  A-u,  a»4.0  km/s 
and  the  wind  speed  is  10  meters/s.  The  results  are  shown  in  figure  4. 

Observe  that  the  linear  relationship  between  R  and  B  persists  throughout  almost 
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FIGURE  4.  THE  THEORETICAL  EFFECTS  OF  "WIND"  NOISE  CONTAMINATION  ON  THE  NOISE 
REDUCING  CAPABILITIES  OF  AN  SI  SYSTEM 
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Che  entire  range  of  B.  Actually  it  is  shown  in  section  5.2  of  the  Appendix 
that  R  will  decrease  as  B  increases  only  for  values  of  1  . 

i+2 

T 

Thus  far,  all  practical  purposes,  equation  48  adequately  described  the  noise 
reducing  capabilities  of  a  vertically  oriented  S-I  system  in  a  field  consisting 
of  Rayleigh  waves  and  wind  noise.  Observe  that  equations  46  and  48  are  iden¬ 
tical.  This  means  in  effect,  that  wind  noise  contamination  is  equivalent  to 
strain  system  noise  contamination  insofar  as  its  impact  on  the  noise  reduction 
capabilities  of  the  vertically  oriented  S-I  system  is  concerned.  This  is  an 
important  result  and  implies  that  the  design  of  the  system  installation  to 
suppress  wind  noise  is  equivalent  in  importance  to  the  design  of  sensors  to 
suppress  system  noise. 

Equation  48  is  not  inherently  limited  to  fields  consisting  only  of  Rayleigh 
wave  noise  and  wind  noise.  It  could  for  example  apply  to  a  field  consisting 
of  two  scattered  Rayleigh  modes  provided  that  Q~^  and  T-*  specify  their  ver¬ 
tical  strain-displacement  response  ratios.  To  illustrate  the  impact  of  model 
interference  we  have  evaluated  equation  47  for  the  cases  where  ^  was 
arbitrarily  assigned  the  value  0.5.  The  results  are  shown  in  figure  5.  The 
important  features  to  note  are  that  R  always  has  a  maximum  somewhere  in  the 
interval  0<B<1  and  vanishes  at  the  end  points  of  this  range.  It  Is  shown  in 
section  5.2  of  the  Appendix  that  the  maximum  occurs  at 

B  =  _ 1_  (49) 

i+a 

T 

and  that  its  value  at  that  point  Is  given  by 


(50) 


Equations  49  and  50  show  that  as  T  I  the  maximum  occurs  at  successively 
smaller  values  of  B  and  diminishes  in  magnitude.  Thus  the  noise  reduction 
capabilities  of  a  strain  inertial  system  can  vary  substantially  depending  upon 
the  strain-displacement  responsibilities  for  the  interfering  modes  and  their 
power  with  respect  to  the  total  power  at  the  output  of  the  strain  sensor. 

2.4.3  Rayleigh  Waves  Contaminated  by  Mantle  P  Waves 

P  waves  arriving  at  near  vertical  angles  of  incidence  are  believed  to  be  a 
significant  constituent  of  the  earth  noise  at  frequencies  greater  than  about 
0.5  to  2.0  Hz.  This  component  is  commonly  referred  to  as  "mantle  P  wave” 
noise.  It  is  important  to  assess  the  effects  that  the  presence  of  a  mantle  P 
wave  component  has  on  the  noise  reducing  capabilities  of  a  vertically  oriented 
S-I  system.  For  our  purposes  it  will  be  convenient  to  express  the  noise 
reduction  R  in  terms  of  the  ratio  of  the  mantle  P  wave  noise  power  to  the 
total  noise  power  as  observed  at  the  output  of  the  inertial  sensor.  If  Bp  is 
this  ratio  then  it  is  shown  in  section  5.3  of  the  Appendix  that  ’ 
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FRACTIONAL  TOWER  OF  CONTAMINATING  MOOE  (OBI 


FIGURE  5.  THE  THEORETICAL  EFFECTS  OF  RAYLEIGH  WAVE  MODE  INTERFERENCE. 
TWO  MOOES  ARE  ASSUMED  TO  BE  PRESENT  AND  9  IS  ASSIGNED  THE 
VALUE  OF  0.5.  T 
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(51) 


('  *(I))bp  «-v 
1  -  (‘ 


where  P  and  V  are  the  moduli  of  the  strain-displacement  response  ratios  for 
the  mantle  P  wave  noise  and  Rayleigh  wave  noise,  respectively.  It  is  also 
shown  in  the  same  section  that  if  the  range  of  velocities  for  the  scattered 
waves  is  small  and  if  the  observation  point  lies  at  surface  of  a  homogeneous 
and  isotropic  half  space  then  (^p2  Is  order  of  10  and 


R  as  B„ 


BP*  1 


Equation  51  has  been  evaluated  explicitly  for  the  case  where  A  ■  p  ,  a»4.0  km/ 
sec  and  the  angles  of  incidence  for  the  incoming  P  waves  lie  in  a  narrow  cone 
centered  on  5°.  The  results  are  shown  in  figure  6.  Observe  that  equation  52 
is  a  good  approximation  to  equation  51  for  most  of  the  range  of  Bp.  It  is 
shown  in  section  5.3  that  R  will  begin  to  deviate  significantly  from  a  simple 
linear  dependence  on  Bp  only  for  values  of  Bp  such  that 

Bp  >  _J_  (53) 

1+P 


at  this  point  R  reaches  a  maximum  value  of 

\l  +H1 

„  _  L  vv/J  (54; 

"max  '  - - p 

1  +  L 
V 

and  thereafter  rapidly  declines  as  Bp  >  1.0.  These  results  may  be  given  the 
following  interpretation.  In  the  interval  0  <  Bp  <  prediction  error 

filtering  will  strip  away  the  Rayleigh  wave  noise  leaving  the  mantle  P  wave 


noise  as  the  residual.  In  the  interval  <_  Bp  <  1.0  prediction  error 

V 

filtering  begin  to  suppress  the  P  wave  noise  and  this  suppression  will  be 
complete  at  Bp  r;  1. 

2.4.4  Scattered  P  Wave  Noise 


In  the  previous  section  it  was  shown  that  when  the  background  noise  consists 
primarily  of  P  waves  arriving  in  a  narrow  cone  surrounding  a  near  vertical 
angle  of  incidence  large  noise  reductions  are  possible.  Let  us  now  consider 
the  limiting  case  where  the  noise  consists  solely  of  P  waves  with  angles  of 
incidence  scattered  in  the  interval  0  <  0  <  where 

0O  -  sin-*n0  (55) 

is  the  maximum  angle  of  incidence  in  the  field.  If  the  strain  and  inertial 
sensors  are  colocated  at  the  surface  of  a  homogeneous  and  isotropic  elastic 
half  space,  then  equations  22  and  26  may  be  used  in  equation  43  to  show  that 
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FIGURES.  THE  THEORETICAL  EFFECTS  OF  "MANTLE  P-WAVE"  NOISE  CONTAMINATION 
ON  THE  NOISE  REDUCING  CAPABILITIES  OF  A  S  I  SYSTEM 


Yei  *  1^1  (n^)|2 

Iee 


Therefore  using  equation  45 


R  -  1  - 


CPI2 


(56) 


(57) 


Equation  57  shows  that  for  the  P  wave  model  under  consideration  the  noise 
reduction  is  an  implicit  function  of  the  maximum  angle  of  incidence  for  the 
scattered  waves.  It  has  been  evaluated  for  representative  values  of  0Q  for 
the  case  where  and  o»4.0  km/s.  The  results  are  shown  in  figure  7. 
Observe  that  in  this  instance,  prediction  error  filtering  will  reduce  the 
noise  by  about  6-9  dB  at  the  most  and  then  only  when  the  maximum  angle  of 
incidence  for  the  waves  in  the  field  is  45°  or  less.  For  larger  maximum 
angles  of  incidence  the  performance  of  the  filtering  process  rapidly  degrades 
to  zero. 


2.5  SIGNAL  DISTORTION 


Let  us  now  consider  the  case  where  the  outputs  of  the  strain  and  inertial 
sensors  consist  of  noise  plus  signal  and  let  us  assume  that  an  optimum  noise 
reducing  filter  has  been  designed  on  the  basis  of  prior  information  about  the 
noise.  It  is  shown  in  section  2  of  the  Appendix  that  the  signal  power  after 
prediction  error  filtering  is  given  by 


S(w) 


hs 

>ii 


1  +  H 


f>ee 

ks 

^ii 


-  2  Re  (H*  <frei) 


KS 

t’ii 


(58) 


3  .3 

where  <t>ee  and  are  the  power  spectral  densities  of  the  signal  at  the  out¬ 
puts  of  the  strain  and  inertial  sensors,  respectively  and  < J>e^  is  the  crosspower 
spectral  density  Jjetween  the  strain  and  inertial  outputs  when  only  the  signal 
is  present.  Let  S(w)  be  the  ratio  of  the  signal  power  after  prediction  error 
processing  to  the  signal  observed  at  the  output  of  the  inertial  sensor.  Then 


S(u>) 


l  +  i  h|  2 


2  Re  (H*  <l>ei) 


s 

ii 


(59) 


S(w)  is  thus  a  measure  of  the  signal  distortion  power  introduced  by  prediction 
error  filtering.  For  the  purposes  of  illustration  let  us  assume  a  plane  P 
wave  signal  and  colocated  sensors  at  the  surface  of  a  homogeneous  and  isotropic 
half  space.  Then  using  the  results  of  section  3  of  the  Appendix  it  can  be 
shown  that 
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FIGURE  7.  THEORETICAL  NOISE  REDUCING 
NOISE  FIELOS  WITH  DIFFERENT 
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(60) 


(61) 


Now  let  us  consider  the  case  of  signal  distortion  caused  by  prediction  error 
filtering  in  a  Rayleigh  wave  noise  field.  It  can  be  seen  from  equation  10 

( ps  _ 

that  in  this  instance  H(oi)  is  real  while  ci  is  a  pure  imaginary.  Therefore 
in  this  case  S(io)  reduces  to  ^ii 


S  (to)  =  1  +  H? 


<t>s 

ee 


(62) 


It  can  be  seen  from  equations  10  and  61  that  S(w)  is  independent  of  frequency 
and  depends  only  on  the  elastic  constants  of  the  medium  and  the  angle  of 
incidence  of  the  signal.  Equation  62  has  been  evaluated  for  the  case  where 
a=4.00  km/s  for  signal  angles  of  incidence  in  the  interval  from  0-90°.  The 
results  are  displayed  in  figure  8.  It-is  important  to  note  that  signal 
distortion  power  is  for  all  practical  purposes  negligible  for  signals  with 
angles  of  incidence  less  than  about  45°.  Since  most  regional  and  teleseismic 
P  wave  signal  satisfy  this  constraint,  the  effects  of  signal  distortion 
arising  from  prediction  error  filtering,  if  the  noise  field  consists  primarily 
of  Rayleigh  waves,  may  be  neglected. 

Now  let  us  consider  the  effects  of  signal  distortion  when  the  ambient  noise 
field  consists  entirely  of  scattered  P  waves.  In  this  instance  it  can  be 
seen  from  equation  28  that  H  Is  a  pure  imaginary.  Hence  equation  59  becomes 


This  function  is  frequency  independent  but  is  sensitive  to  the  range  of  inci¬ 
dent  angles  characterizing  the  scattered  P  wave  noise  field  through  H  and  to 

I 

the  incident  angle  of  the  signal  through  the  ratio  C1  ,  Equation  63  has 

<P  ? 

l  i 

been  evaluated  for  the  ranges  of  noise  field  incident  angles  used  in  section 
2.4.4  and  for  representative  signal  incident  angles  in  the  interval  0  <  0S  <  90°. 
The  results  for  the  case  where  and  u**4.0  km/s  are  shown  in  figure  9.  It 
is  important  to  observe  that  signal  distortion  can  act  to  both  increase  and 
decrease  the  signal  power  observed  at  the  output  of  the  prediction  error 
filter  relative  to  the  power  observed  at  the  output  of  the  inertial  sensor. 
Generally  speaking  these  curves  show  that  if  the  angle  of  incidence  for  the 
signal  lies  outside  the  range  of  incident  angles  for  the  scattered  P  wave  noise 
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FIGURE  8.  THEORETICAL  SIGNAL  DISTORTION  EFFECTS  FOR  RAYLEIGH-WAVE  NOISE 
PREDICTION 
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FIGURE  9.  THEORETICAL  SIGNAL  DISTORTION  EFFECTS  FOR  P  WAVE 
NOISE  PREDICTION 
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the  signal  power  will  be  Increased  by  prediction  error  filtering.  This 
increase  is  substantial  in  those  cases  where  the  angles  of  incidence  is  less 
than  about  12°  and  the  signal  angle  of  Incidence  is  greater  than  about  15°. 
This  result  suggests  that  if  there  are  sites  where  so  called  mantle  P  waves 
are  the  sole  contribution  to  the  noise  field,  the  signal  distortion  effect 
could  provide  significant  Increases  in  regional  P  wave  SNRs.  As  one  might 
expect,  when  the  signal  angle  of  incidence  lies  within  the  range  of  those 
characterizing  the  noise,  the  signal  power  at  the  output  of  the  prediction 
error  filter  will  generally  be  less  than  or  equal  to  the  power  observed  at 
the  output  of  the  inertial  sensor. 


2.6  SNR  ENHANCEMENT 

Let  us  define  the  SNR  enhancement,  I,  as  the  ratio  of  the  power  SNR  observed 
at  the  output  of  the  prediction  error  filter  to  the  power  SNR  observed  at  the 
output  of  the  Inertial  sensor.  It  then  follows  from  the  definitions  of  the 
noise  reduction,  R,  and  the  signal  distortion  S,  that 

I  -  S  (6A) 

R 


Let  us  first  consider  the  case  of  Rayleigh  wave  noise  contaminated  by  strain 
system  noise.  Utilizing  the  results  presented  in  sections  2.4  and  2.5  we 
obtain  the  curve  shown  in  figure  10  where  SNR  enhancement  is  plotted  as  a 
unction  of  signal  angle  of  incidence.  In  this  particular  example,  B,  the 
ratio  of  system  noise  power  to  total  noise  power  at  the  output  of  the  strain 
sensor,  was  chosen  to  be  0.01.  We  know  from  our  previous  discussion  that  the 
effects  of  signal  distortion  may  be  safely  neglected  for  most  regional  and 
teleseismic  P  waves  and  that  R  =  B.  Therefore,  for  all  practical  purposes 

I  =  l  (65) 

B 


In  other  words,  P  wave  signal  enhancement  in  a  Rayleigh  wave  noise  field  is 
totally  controlled  by  the  strain  system  noise  power. 

We  also  know  from  the  results  of  section  2.4.2  that  if  the  Rayleigh  wave 
field  is  contaminated  by  wind  noise  then  R  B,  where,  in  this  instance  B  is 
the  ratio  of  wind  noise  power  to  total  noise  power  observed  at  the  output  of 
the  strain  sensor.  Thus  the  curve  shown  in  figure  10  is  equally  representa¬ 
tion  of  the  SNR  enhancement  obtained  in  a  noise  field  whose  power  is  99% 
Rayleigh  waves  and  1%  wind  noise.  Furthermore,  it  is  clear  that  equation  65 
also  approximately  describes  the  SNR  enhancement  that  prediction  error  filter¬ 
ing  will  provide  in  these  types  of  field.  More  generally  it  can  be  inferred 
from  this  discussion  that  if  Bs  and  Bw  are  the  ratios  of  the  system  noise  and 
wind  noise  powers  to  the  total  noise  observed  at  the  output  of  the  strain 
sensor,  then  the  SNR  enhancement  can  be  written  approximately  as 

I  =  1  (66) 

b8+bw 
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SNR  IMPROVEMENT  IN  dB 
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FIGURE  10.  THEORETICAL  P-WAVE  SNR  ENHANCEMENT  IN  A  RAYLEIGH-WAVE  NOISE  FIELD 
CONTAMINATED  BY  1%  STRAIN  SYSTEM  NOISE  POWER 

G 11260 

-25- 


TR  80-6 


Therefore  Che  SNR  enhancement  is  controlled  by  the  sum  of  the  system  and  wind 
noise  powers  observed  at  the  output  of  the  strain  sensor. 

We  also  observed  in  section  2.4.3  that  if  the  noise  consists  of  scattered 
Rayleigh  waves  contaminated  by  "mantle  P  waves"  then  prediction  error 
filtering  will  reject  the  Rayleigh  wave  noise  leaving  the  mantle  P  wave  com¬ 
ponent  untouched.  In  other  words  the  noise  prediction  filter  predicts  only 
the  Rayleigh  wave  component.  Thus,  based  upon  our  previous  discussion,  the 
effects  of  signal  distortion  may  be  neglected.  Under  these  circumstances  the 
SNR  enhancement  becomes 


1  ;  !_ 

Bp  (67) 

where,  in  this  instance,  B  is  the  ratio  of  the  mantle  P  wave  noise  power  to 
the  total  noise  power  as  observed  at  the  output  of  the  inertial  sensor. 

As  a  final  example  let  us  consider  SNR  enhancement  in  a  noise  field  consisting 
of  scattered  P  waves.  Combining  the  results  shown  in  figures  7  and  9  accord¬ 
ing  to  equation  64  yields  the  SNR  enhancement  data  shown  in  figure  11. 
Comparison  of  these  results  with  the  signal  distortion  curves  shown  in  figure 
9  illustrate  that  for  this  particular  type  of  noise  field  the  SNR  enhancement 
produced  by  prediction  error  filtering  is  significantly  influenced  by  signal 
distortion.  It  is  also  worth  noting  that  the  P  wave  SNR  enhancement  can  be 
substantial  if  the  angle  of  incidence  of  the  signal  is  greater  than  the  angles 
of  incidence  characterizing  the  noise  field. 
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P-WAVE  NOISE  FIELD 


FIGURE  1 1.  THEORETICAL  P-WAVE  SNR  ENHANCEMENT  IN  A  SCATTERED 
P-WAVE  NOISE  FIELD 
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3.  GENERAL  DESIGN  AND  INSTALLATION  CRITERIA 


3.1  SYSTEM  NOISE  CRITERIA 

The  arguments  developed  in  the  previous  section  demonstrate  that  a  vertically 
oriented  S-I  system  has  the  potential  for  wide  band  isotropic  P  wave  SNR 
enhancement.  Furthermore,  they  provide  theoretical  confirmation  of  past 
experimental  observations  which  have  indicated  that  relatively  high  system 
and  wind  noise  levels  at  the  output  of  the  strain  sensor,  as  well  as  com- 
positionally  complex  earth  noise  fields,  have  been  the  principal  factors 
inhibiting  the  realization  of  this  potential.  More  importantly,  however, 
they  provide  the  technical  insight  necessary  to  develop  a  systematic  program 
to  minimize  the  impact  of  these  factors  on  the  P  wave  SNR  enhancement  capabil 
ities  of  future  S-I  systems.  For  example,  as  shown  by  equation  64,  given 
optimum  earth  noise  conditions  (i.e.  scattered  Rayleigh  waves  of  one  mode) 
the  P  wave  SNR  enhancement  capabilities  of  a  vertically  oriented  S-I  system 
will  vary  inversely  in  proportion  to  the  ratio  of  the  system  noise  power  to 
the  total  noise  power  at  the  output  of  the  strain  sensor.  In  other  words,  a 
first  prerequisite  to  the  effective  utilization  of  an  S-I  SNR  enhancement 
system  is  to  maintain  strain  system  noise  levels  low  in  comparison  to  earth 
noise  levels  in  the  bandwidth  of  interest.  This  is  not  a  trivial  problem 
when  one  is  interested  in  the  enhancement  of  P  wave  SNRs  in  the  0.5  to  20  Hz 
bandwidth.  A  few  simple  calculations  are  sufficient  to  illustrate  this  point 
It  is  well  known  that  in  seismically  quiet  locations  that  the  vertical 
displacement  background  noise  level  is  of  the  order  of  a  few  millimicrons  or 
less  at  frequencies  of  1.0  Hz  and  greater.  For  the  purposes  of  demonstration 
let  us  assume  that  this  noise  consists  solely  of  one  Rayleigh  mode  with  a 
phase  velocity  of  1  km/sec.  Then  utilizing  equation  3.34  in  section  3.2  of 
the  Appendix  it  can  be  shown  that  the  corresponding  vertical  strain  noise  in 

the  vicinity  of  one  Hz  will  be  of  the  order  of  10~'*2.  Now  typical  system 

noise  levels  for  vertical  strain  sensors  used  in  previous  experiments  were 
optimistically  estimated  to  be  about  2.5  x  10-^2  in  the  vicinity  of  one  Hz 
(Shopland  and  Kirklin,  1969).  Thus  according  to  equation  64,  the  maximum  P 
wave  SNR  enhancement  that  could  have  been  achieved  under  these  conditions 
could  not  be  more  than  a  few  dB.  An  important  implication  of  this  simple 
illustration  is  that  because  of  the  system  noise  problem  we  have  never  been 
able  to  test  the  S-I  concepts  in  seismically  quiet  regions.  While  it  may  be 

argued,  with  some  justification,  that  the  presumably  complex  structure  of  the 

ambient  earth  noise  in  such  regions  rules  out  the  effective  utilization  of  an 
S-I  system,  the  experimental  data  to  support  this  position  is  almost  totally 
lacking.  Furthermore,  it  is  worth  noting  that  if  the  earth  noise  in  such 
regions  does  consist  solely  of  P  waves  with  near  vertical  angles  of  incidence 
as  some  observers  have  suggested,  then,  as  shown  in  section  2.6,  our  theoret¬ 
ical  calculations  indicate  that  large  regional  P  wave  SNR  gains  may  be 
possible.  In  order  to  experimentally  test  this  possibility  we  must  have  a 
strain  sensor  with  a  system  noise  level  of  10”^  or  less  in  the  0.5  to  20  Hz 
bandwidth. 


3.2  INSTALLATION  DEPTH  CRITERIA 


The  fabrication  of  a  high  sensitivity  strain  sensor  is  only  the  first  step  in 
the  implementation  of  an  effective  S-I  P  wave  SNR  enhancement  system.  It  can 
be  seen  from  the  arguments  presented  in  section  2. A. 2  that  at  a  surface  loca¬ 
tion  wind  noise  and  strain  system  noise  have  essentially  the  same  impact  on 
the  P  wave  SNR  enhancement  capabilities  of  a  vertically  oriented  S-I  system. 
Unfortunately,  as  indicated  by  the  following  argument,  the  wind  noise  problem 
can  be  potentially  far  more  severe  than  the  system  noise  problem.  It  has 
been  shown  both  theoretically  and  experimentally  that  the  earth's  response  to 
the  turbulent  pressure  fluctuations  associated  with  the  *ocal  wind  can  be  a 
major  source  of  seismic  noise.  Consider,  for  example,  an  observation  point 
at  the  surface  of  a  homogeneous  and  isotropic  half  space  for  which  A**P  and  u 
is  of  the  order  of  10^  microbars.  Then  from  equation  5.11  in  section  5.2  of 
the  Appendix  it  can  be  seen  that  a  pressure  fluctuation  of  only  one  microbar 
will  produce  a  vertical  strain  of  about  10"^.  This  is  roughly  2  orders  of 
magnitude  greater  than  the  vertical  strains  associated  with  the  Rayleigh  wave 
noise  discussed  earlier.  Thus,  at  seismically  quiet  sites  where  the  strain 
sensor  is  located  at  or  very  near  the  surface  of  the  earth  the  strain  output 
at  frequencies  greater  than  1.0  Hz  will  be  dominated  by  wind  noise  during 
intervals  of  local  atmospheric  turbulence.  This  phenomenon  was  observed 
during  early  S-I  experiments  at  WMSO  and  while  the  specific  nature  of  the 
problem  was  never  recognized,  the  correct  solution  was  implemented.  They 
chose  to  install  the  vertical  strain  sensor  at  a  depth  of  about  18  meters. 
Utilizing  the  results  obtained  later  by  Sorrells  (1971)  it  can  be  shown  that 
the  attenuation  of  pressure  generated  vertical  strains  with  depth  in  a  half 
space  in  which  A-u  is  given  by 


* 


(  d\  «  (l  +  2  -rr  d 

W  l2  7 


)  exp 


(68) 


where  d  is  the  depth,  £  is  the  convective  wave  length  for  the  pressure  change 
and  tp  is  the  ratio  of  the  vertical  strain  at  depth  d  to  the  vertical  strain 
at  the  surface.  Now  at  frequencies  greater  than  0.5  Hz  the  convective  wave 
lengths  for  wind  generated  pressure  changes  will  probably  rarely  be  greater 
than  20  meters  or  so.  Thus  for  the  WMSO  experiment  fi  would  generally  be  less 
than  or  approximately  equal  to  1.  Therefore 


p  >  +  27tJ  exp  (-277)  (69) 

~  38  dB 


In  other  words,  installation  of  the  vertical  strain  sensor  at  a  depth  of  18 
meters  should  reduce  the  wind  noise  by  roughly  2  orders  of  magnitude  or  more. 
Therefore,  the  wind  noise  at  the  input  to  the  buried  sensor  would  probably  be 
no  more  than  10”12  at  frequencies  greater  0.5  Hz.  This  is  roughly  the  same 
order  of  magnitude  as  the  strain  system  noise.  Thus,  for  all  practical  pur¬ 
poses,  the  shallow  burial  technique  Implemented  during  the  WMSO  experiments 
was  sufficient  to  eliminate  the  wind  noise  from  the  visual  records.  Earlier, 
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we  determined  that  strain  system  noise  levels  on  the  order  of  10“^  in  the 
0. 5-20  Hz  bandwidth  were  needed  in  order  to  test  the  S-I  method  in  seismically 
quiet  regions.  Since  wind  noise  imposes  an  identical  limit  on  the  P  wave  SNR 
enhancement  capabilities  of  a  vertically  oriented  S-I  system  we  must  impose 
the  same  constraint  on  the  maximum  acceptable  wind  noise  level  in  the  band¬ 
width  of  interest.  In  order  to  realize  this  constraint  we  must  install  the 
S-I  system  at  depths  greater  than  the  18  meters  used  in  the  past.  Based  upon 
calculations  utilizing  equation  68  the  minimum  "safe"  depth  for  the  installa¬ 
tion  of  future  systems  should  probably  be  no  less  than  30-40  meters. 


3.3  SYSTEM  DIMENSIONS 

It  should  be  recalled  that  in  our  theoretical  development,  colocated  sensors 
and  point  strain  observations  were  invariably  assumed.  The  practical  implica¬ 
tions  of  these  comments  merits  further  comment.  It  is  a  well  known  experi¬ 
mental  fact  that  if  the  earth  noise  field  is  sampled  at  separate  points  the 
degree  of  coherence  between  the  observations  trends  to  diminish  monotonically 
as  the  separation  distance  Increases.  Thus  as  a  general  rule 

Yit  (*r,co)  <1  6 r  /  0  (70) 


where  6r  is  the  separation  distance  between  the  two  observation  points  and 
the  repeated  index  ii  mean  that  Identical  inertial  component  of  earth  move¬ 
ment  are  sampled  at  the  two  observation  points.  It  can  also  be  shown  that 
for  earth  noise  fields  which  vary  spatially  as  well  as  temporally. 

Yet  (6 r,w)  <  Y^  (6r,w)  8r  /  0  ( 71 ) 


where  the  index,  ei,  denotes  a  strain  component  of  earth  movement  is  observed 
at  one  point  while  an  inertial  component  is  observed  at  the  other.  The 
equality  holds  for  the  special  case  where  the  field  is  a  plane  wave.  Thus 
the  square  of  the  simple  coherence  between  inertial  outputs  separated  by  a 
distance,  *r,  may  be  taken  as  the  upper  bound  on  the  square  of  the  simple 
coherence  between  strain  and  inertial  outputs  separated  by  the  same  distance. 
Now  for  the  purposes  of  illustration,  let  us  assume  that  the  earth  noise 
consists  of  waves  whose  apparent  horizontal  velocities  are  contained  in  the 
interval  «>  <_  c  <_  c  %  and  that  the  power  spectral  densities  of  the  individual 

waves  are  independent  of  their  velocities.  In  this  instance,  it  can  be  shown 
that  the  square  of  the  simple  coherence  between  two  vertical  inertial  outputs 
is  given  by 


Yji  (<5r,w) 


J 


p  u>6r 

1  (~TT> 


oi>^6r^ 


(72) 


where  Jj(£)  is  a  Bessel  function  of  the  first  kind  and  order  1.  Now  it  has 
been  shown  in  section  2.4  that  neglecting  the  effects  of  signal  distortion 
the  SNR  enhancement  capabilities  of  a  vertically  oriented  S-I  system  is  given 
by 
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Therefore  for  the  problem  under  consideration 

I  <  _ l _  (73) 

1  -  **  A  (6 

a)2  6r2 

Now  suppose  we  require  that  the  limitations  on  P  wave  SNR  enhancement  arising 
from  lateral  separation  of  the  strain  and  inertial  sensors  be  no  greater  than 
the  strain  system  noise  and  wind  noise  constraints;  i.e. 

lmLn  1  100  (74) 

Applying  this  condition  yields 

<Sr  <  0.06  (75) 

f 

where  f  is  the  frequency  in  Hz.  Now  there  are  regions  where  Cj,  may  be  of  the 
order  of  a  few  hundred  meters  per  sec  or  less  whereas  we  are  interested  in  P 
wave  SNR  enhancement  at  frequencies  as  high  as  20  Hz.  Thus  in  order  to  avoid 
significant  limitations  on  the  SNR  enhancement  capability  of  an  S-I  system  at 
high  frequencies  which  arise  from  noise  field  decorrelation,  the  separation 
distance  between  the  strain  and  inertial  sensors  should  be  no  more  than  a  few 
meters  at  the  most.  Hence,  strict  adherence  to  the  colocation  assumption  is 
critical  to  the  successful  utilization  of  an  S-I  system  for  regional  P  wave 
SNR  enhancement. 

By  a  similar  line  of  reasoning  it  can  be  shown  that  similar  constraints  apply 
to  the  separation  distance  between  the  sensing  elements  of  the  strainmeter; 
i.e.  if  one  wishes  to  effectively  utilize  the  S-I  method  of  P  wave  SNR  enhan¬ 
cement  at  frequencies  as  high  as  20  Hz  in  noise  fields  where  the  horizontal 
velocities  may  be  as  low  as  a  few  hundred  meters/sec  then  the  strain  and 
length  should  be  no  more  than  a  few  meters  at  the  most.  In  other  words  the 
strain  measurement  should  closely  approximate  a  point  observation. 


3.4  INSTALLATION  NOISE 

At  this  point  in  our  discussion  we  should  logically  consider  the  magnitude 
and  the  means  for  controlling  "installation  noise".  By  this  term  we  mean  that 
which  might  arise  from  sources  Immediately  adjacent  to  the  S-I  system. 
Unfortunately,  we  have  neither  the  experimental  data  nor  the  theory  to  attack 
this  problem  at  this  time  or  even  to  determine  whether  it  will  substantially 
limit  the  P  wave  SNR  enhancement  capabilities  of  an  S-I  system.  Let  us  state 
therefore  that  our  goal  is  to  maintain  strain  installation  noise  at  or  below 
the  10”13  level  in  the  bandwidth  0.5  to  20  Hz.  If  preliminary  experiments 
with  prototype  models  of  future  S-I  systems  demonstrate  that  this  goal  is  not 
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achievable  then  investigations  must  be  carried  out  to  develop  methods  for 
minimizing  the  impact  of  this  type  of  noise  on  the  P  wave  SNR  enhancement 
capabilities  of  an  S-I  system. 


3.5  SUMMARY  OF  DESIGN  AND  INSTALLATION  CRITERIA 

The  principal  design  and  installation  constraints  developed  in  the  previous 
sections  may  be  briefly  summarized  as  follows. 

1.  Strain  System  Noise  Level:  10“^  in  the  0.5-20  Hz  bandwidth. 

2.  Maximum  Linear  Dimensions:  Less  than  3  meters. 

3.  Minimum  Installation  Depth:  40  meters. 

4.  "Installation"  Noise  Level:  10“^  t^e  0.5  -  20  Hz  bandwidths. 

4.  CONCEPTUAL  DESIGN  FOR  AN  ADVANCED  STRAIN- INERTIAL  SEISMOGRAPH  SYSTEM 


4.1  GENERAL  DESCRIPTION 

A  conceptual  design  for  a  vertically  oriented  S-I  system,  capable  of 
installation  in  a  shallow-borehole,  and  which  generally  satisfies  the  design 
constraints  discussed  in  the  previous  section,  has  been  developed.  The  unit, 
which  has  been  designated,  the  Strain-Inertial  Seismograph  System,  Model  52000, 
will  contain  components  to  sense  vertical  displacement  and  differential  ver¬ 
tical  displacement  (vertical  strain)  to  filter,  and  combine  the  strain  and 
inertial  outputs  and  to  display  both  the  individual  and  combined  outputs  as  a 
function  of  time. 

As  shown  by  the  block  diagram  in  figure  12,  4  major  interfaces  conceptually 
divide  the  system  into  3  major  components,  the  borehole,  the  S-I  seismometer 
and  the  wellhead  terminal.  These  components  are  discussed  in  greater  detail 
in  the  following  paragraphs. 


4.2  S-I  BOREHOLE 

The  S-I  borehole  will  be  engineered  to  satisfy  certain  special  conditions. 

In  particular,  as  shown  in  figure  13,  one  section  of  the  borehole  casing  will 
contain  a  precision  instrument  support.  This  support  will  be  welded  to  a 
casing  section  then  stress  relieved  prior  to  entry  into  the  borehole.  It 
will  then  be  joined  to  standard  casing  to  place  the  support  at  a  preselected 
depth. 
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Since  the  support  prohibits  float  shoe/squeegee  cementing  treatment,  the 
cement  will  be  applied  externally  to  the  casing.  The  cement  will  be  chosen 
to  optimize  the  interface  between  the  casing  and  the  local  formation. 

The  minimum  borehole  depth  and  maximum  deviation  angle  will  be  appropriate 
for  high  frequency  P  wave  enhancement  (i.e.,  frequencies  greater  than  1.0  Hz). 

4.3  THE  S-I  SEISMOMETER 

The  S-I  seismometer  consists  of  an  element  which  senses  linear  differential 
displacement  between  the  upper  and  lower  instrument  supports  (the  strain 
seismometer)  and  an  element  which  senses  acceleration  at  the  upper  seat  of 
the  instrument  support  (the  inertial  seismometer).  It  also  contains  elec¬ 
tronics  for  power  conditioning,  communication  control  and  monitoring. 

The  strain  seismometer  consists  of  the  strain  rod  assembly,  the  strain  trans¬ 
ducer  assembly  and  the  associated  electronic  circuitry. 

The  strain  rod  will  be  matched  to  the  borehole  casing  support  and  will 
operate  under  essentially  constant  stress  conditions.  The  strain  transducer 
will  achieve  a  large  linear  travel  range  (+2  mm  typical)  by  using  a  variable¬ 
coupling  type  capacitor  construction.  Desired  resolution  will  be  achieved  by 
operating  4  independent  capacitance  bridges  with  summation  after  amplification. 
Consideration  will  also  be  given  to  organizing  the  bridges  in  quadrants  to 
provide  the  capability  to  sense  horizontal  differential  displacements. 

The  strain  seismometer  electronic  circuitry  will  perform  the  following 
functions : 

a.  Capacitance  bridge  excitation; 

b.  Capacitance  bridge  output  amplification; 

c.  Bridge  signal  summation,  amplification,  and  demodulation; 

d.  Develop  low-gain  dc-coupled  output; 

e.  Develop  high-gain  short  period  output. 

The  high-gain,  short  period  output  will  be  the  primary  signal  output.  Its 

characteristics  will  not  be  altered  by  the  presence  of  large  (10“'  m/m)  long 

term  strain  components.  Depending  on  whether  "quadrant"  or  "linear-only" 
organization  of  the  transducer  assembly  is  utilized,  electronic  circuits  will 
be  provided  for  each  primary  output.  Overall,  the  design  of  the  strain 
seismometer  is  expected  to  minimize  the  effects  of  temperature  gradients, 
pressure  gradients  residual  stresses,  spurious  modes  and  cross  axis  disturb¬ 
ances. 

The  conceptual  design  includes  a  vertical  KS-36000  as  the  inertial  sensor. 

It  will  be  used  without  leveling  mechanism  to  maintain  alignment  with  the 
strain  sensor.  It  can  accommodate  tilts  of  3-4  degrees  at  temperate  latitudes 


and  temperatures.  Special  adjustments  will  have  to  be  made  to  accommodate 
extreme  conditions.  The  inertial  electronic  circuits  will  be  similar  or  com¬ 
mon  with  those  of  the  strain  seismometer. 


It  is  possible  that  two  other  inertial  seismometers,  the  S-700  and  the  44000, 
could  replace  the  KS-36000  in  the  S-I  system.  In  general,  for  these  units  to 
qualify  for  use  in  the  S-I  system,  the  response  should  not  change  over  the 
short  period  band  (0.5-10  Hz)  and  the  short  period  seismometer  noise  should 
not  exceed  3  x  10“^  (m/sec^)^/Hz.  While  the  20171A  (or  23900)  satisfies  the 
noise  constraint,  its  response  and  response  stability  characteristics  make  it 
less  desirable  for  the  S-I  system. 


4.4  THE  WELLHEAD  TERMINAL 

The  S-l  wellhead  terminal  will  perform  the  following  functions 

a.  Power  conditioning 

b.  Operating  controls 

c.  Signal  conditioning 

d.  Data  conversion 

e.  Data  storage  and  retrieval 

f.  Data  processing 

g.  Data  display 


Power  conditioning  will  be  provided  by  the  wellhead  terminal  power  assembly 
which  interfaces  with  the  primary  power  source  at  the  site.  It  will  provide 
power  to  the  terminal  and  the  S-I  seismometer.  A  signal  conditioner  will 
interface  with  the  S-I  seismometer  signal  outputs  and  will  provide  analog 
response  and  shaping  and  anti-alias  filtering  for  primary  data  channels.  The 
conditioned  signals  will  be  applied  to  a  data  converter  and  may  be  selected 
for  display. 

Initially,  the  individual  strain  and  inertial  conditioned  signals  will  be 
digitized  and  recorded  on  magnetic  tape  at  the  site.  The  strain  and  inertial 
data  will  then  be  examined  in  conjunction  with  supplementary  noise  field  data 
to  determine  optimum  P  wave  enhancement  processing  techniques.  After  proces¬ 
sing  definition  specifications  will  be  prepared  for  the  S-I  converter  and 
processor. 
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4.5  SUMMARY  OF  PRIMARY  SYSTEM  CHARACTERISTICS 


The  primary  characteristics  for  the  major  components  of  the  S-I  system  are 
summarized  below: 


Borehole 


Depth 

Casing  diameter 

Slant  angle 

Strain-inertial  seismometer 

Noise  power,  eq,  input 
Strain 
Inertial 

Primary  output  response 
Strain 
Inertial 

Secondary  output 
Strain 
Inertial 

Strain  rod  length 
Wellhead  terminal 


40  meters,  (131  ft) 

7  inch  o.d.,  20  lb/ft,  API 
3-degrees 

4  x  10“26  (m/m)^/Hz 

1  x  10-J8  (m/s^)^/Hz 

1  x  10®  volt/(m/m),  0.2  to  10  Hz 
3  x  10^  volt/(m/s),  0.2  to  10  Hz 

5  x  10^  volt/ (m/m) ,  dc  to  10  Hz 

1  x  10-*  volt/(m/s^),  dc  to  10  Hz 
1  meter 


Characteristics  of  the  wellhead  terminal  assemblies  will  evolve  with  the 
acquisition  and  operating  systems.  Initially,  the  SDCS  acquisition  equip¬ 
ment  can  be  used  with  special  signal  conditioning  circuits. 
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5.  GENERAL  UTILIZATION  CONSTRAINTS 


Development  of  the  S-I  seismograph  system,  whose  conceptual  design  was 
discussed  in  the  previous  section,  will  eliminate  many  of  the  problems  that 
plagued  earlier  attempts  to  reduce  the  S-I  method  of  P  wave  SNR  enhancement 
to  practice.  The  remaining  questions  surrounding  its  ultimate  value  to  the 
underground  nuclear  test  detection  program  relate  to  the  nature  of  the  propa¬ 
gating  earth  noise  field  in  the  0.5-20  Hz.  In  this  regard,  certain  general 
constraints  upon  the  effective  utilization  of  an  S-I  system  are  worth  noting. 
The  S-I  method  was  originally  conceived  to  enhance  P  wave  SNRs  in  the  presence 
of  earth  noise  consisting  of  scattered  Rayleigh  waves  of  a  single  mode.  As 
we  have  seen  from  the  discussion  in  section  2,  it  will  perform  remarkably 
well  under  these  conditions,  being  limited  only  by  the  fractional  power  of 
the  system  and  wind  noise  power  at  the  output  of  the  strain  sensor  in  the 
bandwidth  of  interest.  As  the  earth  noise  field  departs  from  this  simple 
structure  the  SNR  enhancement  capabilities  of  the  S-I  system  will  diminish. 

The  simple  example  of  a  field  consisting  of  Rayleigh  waves  and  "mantle"  P 
waves  succinctly  illustrates  this  point.  We  see  from  the  results  of  section 
2  that  in  this  instance  the  S-I  system  will  simply  strip  away  the  Rayleigh 
waves  leaving  the  mantle  P  wave  noise  intact.  More  generally  we  may  infer 
from  this  example,  that  if  the  noise  field  contains  well  organized  high  and 
low  velocity  components  in  the  bandwidth  of  interest  then  the  fractional 
power  characterizing  the  high  velocity  component  establishes  the  limit  on  the 
P  wave  SNR  capabilities  of  the  S-I  system.  In  other  words  if  the  noise  power 
is  dominated  by  a  high  velocity  component  then  the  P  wave  SNR  enhancement 
capabilities  of  an  S-I  system  will  be  negligible.  The  term  'is  dominated' 
requires  some  further  clarification.  We  have  seen  in  section  2  that  if  the 
earth  noise  consists  entirely  of  P  waves  with  near  vertical  angles  of  inci¬ 
dence,  large  P  wave  SNR  gains  for  regional  signals  are  theoretically  possible. 
Thus  there  exists  a  narrow  power  range  in  which  the  SNR  enhancement  capabili¬ 
ties  may  actually  increase  as  the  percentage  of  high  velocity  earth  noise 
increases.  From  the  results  contained  in  section  2  this  range  is  approximately 
defined  by 


0.999  i  Bp  1  1 


(76) 


where  Bp  is  the  ratio  of  the  high  velocity  noise  power  to  the  total  noise 
power  observed  at  the  output  of  the  inertial  sensor.  Therefore,  the  SNR 
enhancement  capabilities  of  an  S-I  system  will  generally  be  degraded  by  the 
addition  of  a  high  velocity  noise  component  to  a  Rayleigh  wave  field  provided 
that  Bp  lies  in  the  interval. 


0  1  Bp  <  0.999 


(77) 


Outside  this  interval  large  regional  P  wave  gains  may  be  possible.  This 
latter  condition  will  hold,  however,  only  if  the  noise  consists  of  P  waves 
with  near  vertical  angles  of  incidence.  If  the  P  wave  noise  field  contains 
components  with  angles  of  incidence  as  large  as  12°-15°,  the  use  of  an  S-I 
system  can  reduce  rather  than  enhance  regional  P  wave  SNRs.  Thus,  effective 
utilization  of  an  S-I  system  in  a  high  velocity  noise  field  requires  a 
special  set  of  conditions  which  may  be  difficult  to  realize  in  practice. 
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Overall,  then,  we  see  that  even  though  development  of  an  advanced  S-l  system 
will  largely  eliminate  the  design  and  installation  problems  encountered  during 
previous  experiments,  its  P  wave  SNR  enhancement  capabilities  will  still  be 
constrained  by  the  composition  and  state  of  organization  of  the  earth  noise 
in  the  bandwidth  of  interest.  The  general  nature  of  the  earth  noise  in  the 
regional  signal  bandwidth  is  briefly  discussed  in  the  following  section. 


6.  EARTH  NOISE  IN  THE  REGIONAL  SIGNAL  BANDWIDTH 


6.1  REVIEW  OF  EXISTING  DATA 

For  our  purposes,  the  regional  signal  bandwidth  is  considered  to  cover  the 
frequency  range  0.5  <  f  <  20  Hz.  During  the  past  two  decades  the  noise  in 
this  bandwidth  has  been  the  subject  of  numerous  investigations.  In  particular, 
the  need  to  develop  reliable  methods  for  the  enhancement  of  teleseismic 
signals  led  to  the  deployment  of  large  scale  horizontal  arrays  and  deep  bore¬ 
hole  seismograph  systems  during  the  mid-1960s.  The  data  provided  by  these 
sensor  systems  formed  the  basis  of  several  experimental  studies  in  the  late 
1960s  and  early  1970s  which  yielded  a  fairly  comprehensive  qualitative  picture 
of  the  earth  noise  at  frequencies  less  than  2  or  3  Hz.  For  our  purposes  it 
is  worth  noting  that  in  the  interval  0. 5-1.0  studies  of  the  data  from  the 
Montana  LASA  (Toksoz  and  Lacoss,  1968,  Haubrich  and  McKamy,  1969,  Lacoss  et 
al,  1969,  Capon  1969)  indicate  that  the  earth  noise  in  this  band  is  predomin¬ 
antly  body  waves  mixed  with  higher  mode  Rayleigh  waves.  Analysis  of  the  data 
from  deep  borehole  seismograph  systems  is  also  consistent  with  this  description 
(Douze,  1967).  This  noise  appears  to  be  derived  principally  from  oceanic 
sources,  although  there  is  some  evidence  that  local  or  regional  sources  may 
make  significant  contributions  during  quiet  intervals  (Iyer  and  Healy,  1972, 
Capon  1973).  The  relative  power  levels  of  the  individual  components  tends  to 
vary  as  a  function  of  frequency  and  as  a  function  of  time.  In  particular, 
high  body  wave  noise  levels  appear  to  be  correlated  with  storms  at  sea 
(Haubrich  and  McKamy,  1969). 

As  shown  in  figure  14  there  is  a  distinct  change  in  the  slope  of  the  ambient 
earth  noise  spectrum  near  1.0  Hz.  This  slope  change  is  thought  to  mark  the 
transition  from  "oceanic"  to  locally  generated  noise.  Much  less  is  known 
about  the  composition  and  structure  of  this  local  noise.  However,  synthesis 
of  the  results  from  investigations  related  to  underground  nuclear  test  detec¬ 
tion  (Frantil  et  al,  1962,  Douze  1964,  1967,  Robertson,  1965,  Sanford  et  al 
1963,  Fix  1972)  and  from  experimental  assessments  of  the  "ground  noise" 
geothermal  exploration  method  (Clacy,  1968,  Whiteford,  1970,  1975;  Douze  and 
Sorrells,  1974;  Goforth  et  al,  1972,  Iyer  and  Hitchcock  1974,  1976;  Liaw  and 
McEvilly,  1979)  yields  a  hypothesis  concerning  the  major  components  of  the 
earth  noise  at  frequencies  greater  than  one  Hz.  The  hypothetical  power 
spectral  densities  of  these  components  are  schematically  illustrated  in 
figure  13  and  are  discussed  in  greater  detail  in  the  following  paragraphs. 

Curve  1  represents  the  power  spectral  density  of  a  residual  component  commonly 
observed  at  sites  that  are  isolated  from  both  cultural  and  naturally  occurring 
sources.  Typical  features  include  one  or  more  peaks  in  the  2-4  Hz  range 
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followed  by  a  more  or  less  featureless  band  In  which  the  spectral  amplitude 
diminishes  as  the  inverse  square  of  the  frequency.  This  particular  estimate 
was  computed  from  data  collected  at  Queen  Creek,  Arizona  (Fix,  1972).  Very 
similar  spectral  envelopes  characterize  the  noise  observed  at  Kanab,  Utah 
(Frantii  et  al,  1962)  and  at  the  Ogdens burg  mine  during  quiet  intervals 
(Isacks  and  Oliver,  1962).  It  has  been  suggested  that  the  residual  component 
consists  primarily  of  body  waves  arriving  at  near  vertical  angles  of  incidence 
(so  called  mantle  P  waves).  However,  this  explanation  is  based  largely  on 
limited  observations  at  frequencies  at  less  than  2-3  Hz.  Thus,  in  our 
opinion,  the  composition  of  this  component  is  largely  undefined  at  this  time. 

At  sites  where  nearby  cultural  or  naturally  occurring  sources  are  intermit¬ 
tently  active  the  residual  component  is  commonly  overridden  by  a  second  com¬ 
ponent  whose  spectral  density  is  schematically  illustrated  by  curve  2  in 
figure  14.  This  component  will  be  referred  to  as  the  surface  wave  component 
because  it  is  believed  to  consist  primarily  of  fundamental  mode  Rayleigh 
waves.  However,  there  is  evidence  to  suggest  that  the  individual  peaks  in 
its  power  spectral  density  may  be  associated  with  the  local  group  velocity 
minima  of  higher  mode  Rayleigh  waves  (Asten,  1978).  The  shape  of  the  spectral 
envelope  seems  to  be  strongly  dependent  on  the  local  velocity  structure  in 
the  uppermost  few  hundred  meters  of  the  earths  crust  and  in  some  instances 
may  be  largely  determined  by  the  structure  in  the  uppermost  few  10s  of  meters 
(Liaw  and  McEvilly,  1979).  The  time  dependence  and  intensity  of  the  noise 
sources  may  also  play  a  significant  role  in  determining  the  local  spectral 
envelope  of  this  second  component.  However,  source  effects  have  been  clearly 
documented  in  only  a  few  isolated  instances  (cf.  Sanford,  et  al,  1968). 

A  third,  wind  related,  component  is  often  observed  at  surface  sites  when  local 
wind  speeds  are  in  excess  of  a  few  miles  per  hour  or  so.  The  power  spectral 
density  is  schematically  illustrated  by  curve  3  in  figure  14.  The  results  of 
experimental  studies  by  Douze  (1964)  coupled  with  studies  by  Sorrells  (1971) 
strongly  suggests  that  this  component  consists  principally  of  near  field, 
quasi-static  deformations  of  the  surface  of  the  earth  in  response  to  wind 
generated  atmospheric  pressure  fluctuations.  The  power  in  this  component  is 
functionally  dependent  upon  the  local  atmospheric  pressure  power  spectral 
density  and  upon  the  local  medium  response  to  pressure  changes  acting  at  the 
surface  of  the  earth. 


6.2  IMPLICATIONS 

The  existing  data  base  concerning  earth  noise  in  the  regional  signal  bandwidth 
is  deficient  for  our  purposes  in  that  it  fails  to  provide  quantitative  infor¬ 
mation  concerning  relative  composition  and  state  of  organization.  As  we  have 
seen  from  our  previous  discussions,  the  P  wave  SNR  enhancement  capabilities 
of  any  future  S-I  system  will  be  critically  dependent  upon  these  factors. 

Thus,  a  complete  assessment  of  the  feasibility  of  developing  an  advanced  S-I 
cannot  be  made  at  ths  time.  Nevertheless,  certain  generalizations  are 
possible.  For  example,  based  upon  the  available  data  it  seems  very  likely 
that  an  S-I  system  can  be  used  to  enhance  P  wave  SNRs  in  regions  where  the 
surface  wave  component  is  predominant.  Significant  SNR  gains  are  likely  to 
occur  only  at  frequencies  greater  than  one  Hz  and  may  vary  substantially 
within  the  1-20  Hz  passband  depending  upon  the  relative  composition  of  the 


-42- 


TR  80-6 


noise  power  at  different  frequencies.  The  relative  composition  of  the  noise 
may  also  vary  substantially  in  time  as  well  as  over  relatively  short  lateral 
distances  depending  upon  the  local  near  surface  geology.  Thus  the  P  wave  SNR 
enhancement  capabilities  of  an  S-I  system  is  very  likely  to  be  strongly  site 
and  time  dependent  in  regions  where  the  surface  wave  component  is  predominant 
in  the  regional  signal  bandwidth.  At  locations  where  the  residual  component 
predominates,  the  P  wave  SNR  enhancement  that  could  be  obtained  through  the 
use  of  an  S-I  system  must  be  termed  highly  speculative  at  this  time.  While 
our  calculations  indicate  that  large  regional  P  wave  SNR  gains  are  possible, 
if  the  residual  noise  does  indeed  consist  of  "mantle"  P  waves,  conclusion 
evidence  to  support  this  description  is  not  currently  available. 

Overall,  the  results  of  the  review  of  the  existing  data  base  demonstrate  that 
the  P  wave  SNR  enhancement  capabilities  may  vary  from  excellent  to  poor 
depending  upon  the  nature  of  the  local  earth  noise  field.  Thus  intelligent 
use  of  future  S-I  systems  will  require  a  priori  information  about  the  local 
properties  of  the  earth  noise  field  In  the  regions  where  they  are  to  be 
deployed. 


6.3  PRELIMINARY  EXPERIMENTAL  RESULTS 

From  the  foregoing  discussion  it  should  be  clear  that  quantitative  information 
concerning  the  composition  and  state  of  organization  of  the  earth  noise  is 
necessary  in  order  to  realistically  assess  the  P  wave  SNR  enhancement  capabil¬ 
ities  of  an  S-I  system.  It  is  equally  clear  that  the  existing  data  base  fails 
to  provide  the  required  levels  of  detail.  Therefore,  as  a  part  of  this  pro¬ 
ject,  preliminary  investigations  were  undertaken  to  expand  the  data  base  per¬ 
taining  to  the  earth  noise  in  the  regional  signal  bandwidth.  During  the 
course  of  these  investigations  we  have  adapted  Liaws  (Liaw,  1978)  method  for 
the  estimation  of  high  resolution  frequency  wave  number  (f-k)  spectra  to  run 
on  the  Vax  11/780  system  and  have  developed  a  working  version  of  Abo  Zena's 
(Abo  Zena,  1979)  method  for  calculating  the  phase  and  group  velocities 
characterising  a  horizontally  layered  elastic  solid.  Abo  Zena's  method  is 
particularly  useful  to  us  in  that  it  provides  us  with  the  capability  to 
accurately  calculate  the  phase  and  group  velocities  for  high  frequency 
Rayleigh  wave  modes  in  thinly  bedded  media.  In  addition,  preliminary  field 
measurements  were  made  at  McKinney,  Texas,  in  order  to  develop  criteria  for 
the  design  of  future  experiments  at  this  site.  The  data  collection  system 
used  to  obtain  the  preliminary  information  was  furnished  by  the  SDCS  program 
sponsored  by  AFTAC/VSC.  A  block  diagram  of  the  system  is  shown  in  figure  15. 
Basically  the  system  consisted  of  six  short  period  seismometers  and  signal 
conditioning  electronics  to  yield  samples  of  the  earth  noise  in  the  0.5-20  Hz 
bandwidth.  The  overall  system  response  to  a  constant  displacement  input  is 
shown  in  figure  16.  In  addition,  a  trigger  mechanism  was  provided  so  that 
data  were  collected  by  the  system  only  when  the  local  wind  speeds  were  less 
than  5  mph.  When  this  condition  was  met,  data  from  the  6  seismographs  were 
sampled  at  125  ms  intervals  and  stored  on  magnetic  tape  for  future  analysis. 
The  sensors  were  deployed  in  a  linear  array.  During  the  measurements  program 
the  interval  between  sensors  was  changed  to  produce  arrays  which  varied  from 
50  to  250  meters  in  length.  The  sensors  were  buried  at  a  shallow-depth  and 
care  was  taken  to  ensure  that  the  local  geology  was  uniform  at  each  site  in 
the  array.  An  example  of  some  of  the  preliminary  results  are  shown  in  figure 
17.  These  are  records  from  the  250  meter  array  during  an  interval  when  the 
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FIGURE  16.  FREQUENCY  RESPONSE  OF  NOISE  DATA  COLLECTION  SYSTEM  TO  A  CONSTANT 
DISPLACEMENT  INPUT 
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FIGURE  17.  REPRESENTATIVE  NOISE  RECORDS  AT  McKINNEY,  TEXAS,  DURING 
A  QUIET  INTERVAL.  INTERVAL  BETWEEN  SEISMOMETERS  IS 
50  METERS. 
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wind  speed  was  less  than  5  raph.  In  addition  the  sample  was  taken  when  the 
effects  of  local  cultural  activity  should  be  at  a  minimum.  Observe  that  the 
records  are  dominated  by  a  quasi-sinusoidal  component  with  frequencies  in  the 
2-4  Hz  range.  This  component  appears  to  be  fairly  well  organized  over  the 
dimensions  of  the  array.  Response  corrected  spectral  estimates  for  the  data 
recorded  at  stations  1  and  6  (the  end  points  of  the  array)  are  compared  for 
reference  to  the  Queen  Creek  spectral  estimate  in  figure  18-a.  The  estimated 
square  of  the  simple  coherence  between  the  data  recorded  at  stations  1  and  6 
is  shown  in  figure  18-b.  Notice  the  coherence  minimum  which  occurs  near  1.5 
Hz.  We  suspect  that  this  marks  the  transition  from  the  body  waves  charac¬ 
terizing  the  oceanic  component  to  the  local  noise  which  we  believe  in  this 
case  to  be  dominated  by  surface  waves.  As  implied  hy  the  time  series,  the 
noise  field  power  is  dominated  by  a  well  organized  2-4  Hz  component.  In 
addition  there  is  some  indication  of  weak,  organized  components  in  the  5-8  Hz 
band  and  near  10  Hz.  The  difference  between  the  McKinney  and  Queen  Creek 
spectra  demonstrates  that  even  under  optimum  environmental  conditions, 

McKinney  is  a  relatively  noisy  site.  In  addition,  we  believe  that  future 
studies  will  reveal  that  a  significant  fraction  of  this  difference  is  caused 
by  the  presence  of  low  velocity  fundamental  and  higher  mode  Rayleigh  waves. 

If  this  is  true,  then  McKinney  would  qualify  as  a  potential  site  for  testing 
the  P  wave  SNR  enhancement  capabilities  of  future  S-l  systems. 

Further  studies  of  the  noise  at  McKinney  should  be  directed  towards  testing 
the  hypothesis  that  a  significant  fraction  of  the  earth  noise  in  the  regional 
signal  bandwidth  consists  of  scattered  fundamental  and  higher  mode  Rayleigh 
waves.  This  will  require  the  deployment  of  two  dimensional  spatial  arrays  of 
sensors  and  the  computation  of  (f-k)  power  spectral  density  estimates.  In 
order  to  assist  in  the  design  of  the  data  collection  array  we  have  made  esti¬ 
mates  of  the  variation  of  coherence  as  a  function  of  separation  distance. 

Some  preliminary  results  are  shown  in  tigure  19.  Notice  that  the  square  of 
the  estimated  coherence  diminishes  more  or  less  monotonical ly  both  as  a  func¬ 
tion  of  increasing  frequency  and  separation  distance.  For  all  practical  pur¬ 
poses  the  field  at  frequencies  greater  than  about  2  Hz  is  seen  to  be  completely 
disorganized  at  separation  intervals  on  the  order  of  a  few  hundred  meters  or 
so.  This  suggests  that  arrays  with  relatively  small  lateral  dimensions  can 
be  used  to  collect  the  data  necessary  for  reliable  (f-k)  spectral  estimates 
of  the  high  frequency  noise  in  the  regional  signal  bandwidth. 
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FIGURE  18.  (a)  REPRESENTATIVE  NOISE  POWER  SPECTRAL  DENSITY  ESTIMATES 
AT  McKINNEY.  TEXAS.  DURING  A  QUIET  INTERVAL.  QUEEN  CREEK 
NOISE  ESTIMATE  .  FIX,  1972  IS  INCLUDED  FOR  REFERENCE. 


frequency  <H«) 

FIGURE  18.  (b)  REPRESENTATIVE  COHERENCE  ESTIMATES  FOR  STATIONS  250 
METERS  APART 
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7.  CONCLUSIONS 


The  results  of  our  Investigation  generally  demonstrate  that  a  vertically 
oriented  S-I  seismograph  system  has  the  potential  for  wide  band,  isotropic  P 
wave  SNR  enhancement.  Failure  to  realize  this  potential  may  be  attributed 
largely  to  deficiencies  in  system  design  and  installation  procedures  which 
placed  severe  limitations  on  the  performance  of  previously  developed  systems. 
These  limitations  can  be  avoided  by  developing  an  advanced  S-I  system  which 
incorporates  the  principal  elements  of  the  conceptual  design  described  in 
this  document.  However,  minimization  of  the  design  and  installation 
constraints  will  not  guarantee  satisfactory  performance  from  an  S-I  system 
for  its  P  wave  SNR  enhancement  capability  is  critically  dependent  on  the  com¬ 
position  and  state  of  organization  of  the  earth  noise  in  the  bandwidth  of 
interest.  Thus,  a  complete  assessment  of  S-I  system  P  wave  SNR  enhancement 
capabilities  requires  quantitative  information  regarding  these  parameters. 

Data  to  the  required  level  of  detail  are  not  generally  available,  particularly 
in  the  regional  signal  bandwidth.  Nevertheless,  it  is  worthwhile  to  speculate 
upon  the  regional  P  wave  SNR  enhancement  capabilities  of  an  S-I  seismograph 
system  based  upon  the  information  available  to  us  at  this  time.  In  particular, 
P  wave  SNR  enhancement  is  unlikely  to  be  impressive  in  the  "oceanic”  noise 
band  (-0. 5-1.0  Hz)  except  perhaps  at  sites  located  near  coastlines,  because 
of  the  apparently  high  percentage  of  P  wave  noise  power.  As  our  calculations 
have  indicated,  the  optimum  S-I  prediction  error  processor  will  simply  strip 
away  the  surface  wave  noise  leaving  the  P  wave  noise  as  a  residual.  In  the 
"local"  noise  band  (-1-20  Hz)  P  wave  SNR  enhancement  Is  likely  to  be  highly 
site  specific  and  can  be  expected  to  vary  substantially  both  as  a  function  of 
time  and  frequency.  At  sites  which  are  remote  from  both  cultural  and  naturally 
occurring  noise  sources  P  wave  SNR  enhancement  is  generally  expected  to  be 
negligible.  However,  the  possibility  that  the  noise  at  these  sites  consist 
of  P  waves  at  near  vertical  angles  of  incidence  (mantle  P  waves)  is  worth 
checking  out,  for  under  these  conditions  large  P  wave  SNRs  are  theoretically 
possible.  On  the  other  hand,  regional  P  wave  SNR  enhancement  may  well  be 
substantial  at  sites  where  nearby  noise  sources  are  active  and  the  local 
geology  will  permit  the  existence  of  high  frequency  Rayleigh  waves.  Under 
these  conditions  the  use  of  an  S-I  system  could  conceivably  transform  a  so- 
called  seismically  noisy  site  into  a  moderately  quiet  one,  at  a  relatively 
low  cost. 

Realization  of  this  capability  would  greatly  increase  the  flexibility  of  the 
deployment  strategy  underlying  the  installation  of  a  regional  underground 
nuclear  test  monitoring  system.  For  this  reason  we  conclude  that  development 
of  an  advanced  S-I  system  is  sufficiently  important  as  to  justify  the  fabri¬ 
cation  of  one  or  more  engineering  models  for  the  purpose  of  additional  evalu¬ 
ation  of  this  particular  P  wave  SNR  enhancement  technique. 
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STRAIN- INERTIAL  SNR  ENHANCEMENT  USING  A  PREDICTION  ERROR  FILTER 


STRAIN- INERTIAL  SNR  ENHANCEMENT  USING  A  PREDICTION  ERROR  FILTER 


1.  FILTER  DESIGN 


Let  n^(t)  and  ne(t)  be  the  outputs  of  colocated  inertial  and  strain  sensors 
when  only  background  noise  is  present.  We  wish  to  apply  a  filter  to  the  strain 
output  in  order  to  predict  the  background  noise  on  the  inertial  transducer. 


Let  ni(t)  be  the  prediction  of  ni(t).  Then 


nj  (t) 


h  (t)  ne  (t-T)  dr 


where  h(-r)  is  the  impulse  response  of  the  prediction  filter.  We  will  choose  the 
filter  properties  such  that 


E  j(  ni  (t)  -ni  (t))2  | 


minimum 


Applying  this  constraint  in  the  frequency  domain  requires  that  the  function 

E  (W)  =  (U)  ♦  |  H  (u>)  | 2  ♦Je  (to)  -  II*  (u)  (w)  -  H  (u)  *  0*0. 

1 . 3 

be  minimized,  where  <f>^  (w)  are  the  power  densities  of  the  noise  at  the 
output  of  the  strain  and  inertial  sensors,  ()>n.  (w)  is  the  cross  power  spectral 
density  relating  the  two  outputs  and  H  (u>)  is1the  frequency  response  of  the 
desired  filter.  Applying  this  constraint  yields 

(u)  *  Ci  (“0  “  H*  (“)  (<*0  -  11  («)  *n  (w)  =  0 

ii  ee  ,  . 

1 .4 

-♦Ji  (<*>)  +  j+JJj*  (u)  -  JH*  (u>)  (u>)  -  jH  (to)  (u>)  =  o, 

which  have  the  solution 


H  (w)  = 


*ej  (u)) 
*ee 
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Substitution  of  eq.  1.5  into  eq.  1.3  yields 


E 


O)  =  4>n.  (to) 
11 


|  4>n  (u>)  |2 

ei 


$n. (w) 

li 


<t>n  (w) 
ee 


l 


=  4>n  (w) 

ii 


{■ 


1.6 


where  y2^  (u)  is  the  square  of  the  simple  coherence  between  the  strain  and 
inertial  outputs.  E  (w)  in  this  case  is  seen  to  be  the  power  spectral  density 
of  noise  remaining  after  subtracting  rij  (t)  from  nj  (t).  Clearly  if  the  noise 
at  the  output  of  the  strain  and  inertial  transducers  is  perfectly  coherent  the 
inertial  noise  can  be  completely  eliminated  by  this  process. 
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2 .  SNR  ENHANCEMENT 


Let  yi  (t)  and  (t)  be  the  outputs  of  the  inertial  and  strain  sensors  when 
both  noise  and  signal  are  present.  Then 


yi  (t)  =  sA  (t)  +  nA  (t) 
ye  (t)  =  se  (t)  +  ni  (t) 

Let  yi  (t)  be  the  strain  output  after  application  of  the  noise  prediction 
filter,  i.e. 


yi 


/•oo 

(t)  =  /  h  (x)  ye  (t-x)  dx 


2.2 


/ 


/ 


h  (x)  se  (t-x)  dx  +  I  h  (x)  ni  (t-x)  dx 


Let  e  (t)  be  the  difference  between  the  inertial  output  and  the  filtered  strain 
output .  Then 


L 


e  (t)  =  Si  (t)  -  /  h  (x)  s0  (t-x)  dx 


-  «*>  ■/: 


h  (x)  ne  (t-x)  dT 


2.3 


Let  4>^ .  (w)  and  (w)  be  the  power  spectra  of  the  signal  and  the  outputs  of 

the  inertial  and  strain  sensors,  respectively.  Let  <f>0^  (u)  be  the  crosspower 
spectrum  between  the  signal  at  the  strain  and  inertial  outputs.  Then  the  power 
spectrum  of  the  signal  after  completion  of  prediction  error  processing  is 


S  (u>)  =  *5.  (u)  +  |ll  M|2  *|e  («o)  -  H*  (u>)  $S.  (u)  .  H  (a,)  *s*  (u) 


<f>S  (w) 
ii 


1  +  I H  (a))  | 2  ^eeM  -  II*  (w)  *ci(a,)  -  H(u)  ^i(a,) 


(“) 

ii 


(«) 


Cl 

<t>? .  (w) 

ii 


2.4 
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From  eq.  1.6  the  power  spectrum  of  the  noise  after  prediction  error  processing 
is 

N (w)  =  E(u>)  =  ♦".(«)  2.5 

If  I (w)  is  the  improvement  in  SNR  after  prediction  error  processing  then 

■  (,)  .  iR.  /  ♦!><“) 

N(w)  j  (ul 

1  +  ^ei(u)  *ie(w)  -  2Re  /  ^ei  ^  *ii  (w)  \ 

*ge(u)  \*ge  (<*>)'  (u)/ 

2.6 
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3.  CONCEPTUAL  MODELS 


3.1  P  WAVE  SIGNALS  IN  A  HOMOGENEOUS,  ISOTROPIC  ELASTIC  HALF  SPACE 

Let  A  and  u  be  the  constants  specifying  the  elastic  properties  of  the  half 
space  and  assume  a  plane  P  wave  is  incident  at  the  free  surface.  Let  t  and 
0  be  the  velocity  potentials  satisfying  the  wave  equations. 


1 

«2C 

<5 t2 

V2i|> 

1 

if 

6 t2 

3. 1 

where  a  and 

8  are  the  P  and  S  velocities 

in  the  medium. 

For 

an  incident  P 

wave  these 

potentials  have  the  form 

r. 

=  exp  [jk  ( ya  z-x 

)]  +  a2 

exp  [jk  (-Ya  : 

exp  (juit)  3.2 

* 

=  ^2  expfjk  ("Tg  : ~ x 

)]  exp 

(jut) 

where  k  is 

the  apparent  horizontal 

wave  number,  z  and  x 

arc 

Cartes  ian 

coordinates 

normal  and  parallel  to 

the  free  surface  and 

=  ( 

(£)2  - 
.  a 

3 . 3 

1/2 

>8  =  ( 

(£)2  - 
■  B 

>) 

3.4 

c  is  the  apparent  horizontal  velocity  which  is  defined  by 


sino 

where  9  is  the  angle  of  incidence. 


3.5 


We  are  principally  interested  in  the  vertical  displacements  and  vertical 
strains.  The  vertical  displacement,  W,  and  the  vertical  strain  are 
related  to  the  potentials  c,  0  through  62 


w 

=  £  ♦ 

6  2 

64/  3.6 

6  x 

6W 

=  $!§.♦ 

3.7 

&  z 

fiz2 

fix2 

while  the  normal  and  tangential 

stresses 

acting  in  planes  parallel  to  the  free 

surface  are  related  to  the  potentials  through 

62c  .  fi20 

0  „  =  u  — -  +  — i-  - 

fi  2  0 

3.8 

-x  fixfiz  fix2 

6z2 

0  =  X  V2j;  +  2p  ^  + 

3.9 

6  z2 

fixfiz 

Imposing  the  condition  that  the  stresses  must  vanish  at  the  free  surface 
results  in  the  following  expressions  for  the  reflection  coefficients  A2  and 
• 


3. 10 


3.11 

3.  12 


3.13 
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Using  equations  3.1,  3.2,  and  3.10-3.13  in  equations  3.6  and  3.7  the 
expressions  for  the  vertical  displacements  and  vertical  strains  are  found 
to  be 


IV  =  jkAjj  Ya 

^exp  (jkyazl  - 

T  exp  (- jkyaz ) 

j  -0  exp  (- jkygz)  |  exp  (j(tot-kx)} 

3.14 

fiW  ...  1  .?  [ 

~  =  k-Ajj  >n  | 

CXP  (jk>azl  + 

P  exp  (-jkyazl 

+  C  exp  (-jkyazl  j  exp^j  (u,t-kxl) 

3.13 

3.2  RAY  1.1: l fJ I  WAVES  IN  A  HOMOGENEOUS,  ISOTROPIC,  ELASTIC  HALE  SPACE 

For  the  purposes  of  this  investigation  it  is  sufficient  to  restrict  our 
attention  to  the  vertical  displacements  and  strains  created  by  surface  point 
loads.  Pol  lowing  the  development  in  Ewing  et  al  (10341  let  us  choose  a 
cylindrical  coordinate  system  (r,  0,  cl  and  displacement  potentials  of  the 
form 


c  =  Ac ' v -  II-;  (krl 
if-  =  Be"v Ho(kr) 


where 


(v'l-  =  k-'-  kr 

i' 


k' 

t> 


10' 

lV 


k- 

3 


3.17 


3.  IS 


and  (krl  is  a  ilankcl  function  of  the  second  kind  and  zero  order. 


Let  q  and  W  be  the  radial 
terms  of  the  potentials 

‘I 

W 


and  vertical  displacements 


3-  V' 

fir 

firfi  z 

fir. 

li  j 

^r  fiif- 

r  fir' 

(  fir 

respectively.  Then  in 


3. 10 
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Now  the  normal  and  tangential  stresses  acting  on  planes  parallel  to  the  free 
surface  are  related  to  the  displacements  through 


where 


3.20 


3.21 


Using  equations  3.16  through  3.19  in  equation  3.20  yields 


o2r  =  u  [2kvA  -  k  (2k2 -k2 g) B]  H2  (kr) 

3.22 

azz  =  M  [ (2k2 -k2)  A  -  2k2v'B]  H2  (kr) , 

at  z  =  0 

Now  suppose  a  vertical  point  force  of  strength  ZH2  (kr)  is  applied  at  the 
origin.  Then 

°zrlz=0  =  0 

3.23 

°zz]  =  Z  H2  (kr) 
z=0  0 

and  equations  3.22  become 

-2  vA  +  2  (k2 -kg)  B  =  0 

^  5.24 

(2k2 -k2 ) A-2k2  v  'B  =  - 
8  p 
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Solving  equations  3.24  leads  to  the  following  expressions  for  A  and  B 


A  = 


(2k2-k2)  z 
F(k)  u 


B  = 


F(k)  u 


where 


Therefore 


F(k)  =  (2k?-kJ)2  .  4k‘,vv ' 


v  (2kz-k-')c'vz  +  2kV  e‘v  z 

B  ’  ■' 

W  =  -  -  l!2(kr) 

F(k)  0 


and 


dw  (2k?-k?)  -  2k:> v 1  v  e“v  =  - 

P  ~"o  (kr) 


dZ 


F(k) 


u 


3.25 


3 . 2b 


3.27 


In  order  to  obtain  the  point  force  solutions  2  is  replaced  by  - 
in  equations  3.27  and  the  resulting  expressions  are  integrated 
with  respect  to  k  in  the  interval  (Kk  <  i.e. 


IV 


dW 

dZ 


l(2k?-k£)  e_v=  ♦  2k-'  e"v'r-l 

TO 


!i'0  (kr)  kdk 


[v  (2k-’-kp  e’vz  -2k-1  v 1  e~  2]  . 

-  0 

F  ( k ) 


(kr)  kdk 


3.28 
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Let  WR  and  ^WR  be  the  Rayleigh  wave  vertical  displacements  and  strains  and 

TT 

let  KR  he  a  real  solution  to  the  equation 

F (k)  =0  3, 


Wi>  =  iL 


Krv  [(2Kr2  -  k2)  e'vz  ♦  2K2  e^'z]  |H2  (KRr) 

( 


3WR  _  -jL 
32  u 


Krv  [u(2K2  -  k2)  e~v~  -2K2  v>  e-v‘z] 


where  now 


(K2  -  k2«)1/2 

K 


"4 '  ki,1/2 


and  F '  (Kr1 


i  "<k,|k 


4K„  1K-’  -  2k2  -  (2v/  +  K2  VJ  +  k2  iL  ) 
R  K  P  R  v  R  v' 


For  the  purposes  of  this  investigation  we  arc  principally  interested  in  t ho 
contributions  from  remote  sources.  Thus  for  KRr  1 


!l;  (KrD 


,  v/~ 


exp  i1!-.  exp  l-j  KRr) 
4 


w«  =  i 


/if  L  'l  U[‘:K"  -  « 

W  r^V  I-UKr) 


,  c-”:  •  2K-  ,  v  :]| 


[-'(k"rT)] 


3:  ~  i'  V*-1 


Kr1/Jv  [vU’kR-kj^)  e~vZ  -  2K2V1  e~v  z]| 


F  (Kp) 
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4.  POWER  AND  CROSS  POWER  SPECTRA  FOR  IDEALIZED  NOISE  FIELDS 


4.1  SCATTERED  P  WAVE  NOISE  FIELD 

Let  us  consider  a  noise  field  that  consists  solely  of  scattered  P  waves  with 
apparent  horizontal  phase  velocities  ranging  from  °°  to  cmin  where 

cmin  5  a  4.1 

Let  us  further  require  that  the  incident  P  wave  field  be  laterally  homogeneous, 
isotropic,  and  that  its  potential  amplitude  be  independent  of  the  apparent 
horizontal  phase  velocity.  If  Nqq  (ic.io.z)  is  the  3  dimensional  displacement 
power  spectrum  of  this  field  then 

Nnn  (*\w,z)=|Aj  CO  ,10)  | ?  |D  (k  ,z1 1  ? 

=  0 


|  k  |  < 

cmin  4.2 

|k|  < 

c  . 
m  1  n 


where  from  cq.  3.14 


D  (k,z)  =  j k  <  Ya  [CXP  (jk\az)  -  T  exp  (-jkyaz)]-  S2  exp  i-jky^z) 


4.3 


and  from  equations  3.10  through  3.13 


4k2  (k2  -  k?)1/2  (k2  -  k2)1/2  -  (k2  -  2k2)2 

a _  g _ _ _ _g _ _ 

4k2  (V  -  k2)l/2  fk2  .  k2 ) 1/2  +  (k2  -  2k2)2 

4k  (k2  -  k2i 1/2  (k2  -  2k-  ) 

_ _  g _ 

4k2  (k2  -  k2)1/2(k2  -k2)1/2  .  (k2  .  2k2)2' 
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Similarly  if  Ngg  (k,u,z)  is  the  3  dimensional  vertical  strain  power  spectrum 
then 


Ngg  (k ,oj,z) 


|Ai(e,oi)  |2  |S  (k,z)  l2  [H  <  “  .  ... 

cmin 


0 


|k|  > 


0) _ 

cmin 


4.6 


where  from  equation  3.15 


6D 

<Sz 


(£,z)=  -  |(k2  -  k2)  ^exp  [j  (k2  -  k2)1/2  Z]  +  r  exp  [-j  Cka  '  k2)1/2z|j 
+  k  (k2  -  k2)1/2^exp  [-j  (k2  -  k2)  z]^  j 


4.7 


Finally  if  NEQ  (£,w, z)  is  the  3  dimensional  vertical  strain-vertical 
displacement  cross  power  spectrum  then 


^(k,z)  D* (k, z) 

1  k 

6z 

=  0 

Ik 

min 


cmin 


4.8 


4.2  RAYLEIGH  WAVE,  NOISE  FROM  SCATTERED,  REMOTE  SURFACE  SOURCES 

Let  r  be  the  radius  vector  from  an  arbitrarily  chosen  origin  to  an  observation 
point.  Let  "r '  be  the  radius  vector  from  the  origin  to  a  point  surface  source 
of  Rayleigh  waves.  From  equation  3.34  the  Fourier  time  transform  of  the  vertical 
displacement  component  of  the  Rayleigh  wave 

WR  (r,r"  w;Kr)  =  p^yi/2  °  (KR’z)  exp  exp  ("jKR  4-9 
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Let  us  now  impose  the  restriction  that  |r'|  >>  |r|.  Then 

Kr  | r-r " |  s  KRr'-KRr  cos  (0-6")  4.10 

and 

1  -  1  4.11 

_ 1 r-r-| !/2  '  (r ") 1/2 

where  r,8  are  the  polar  coordinates  of  the  observation  point  and  r',0'  are 
the  polar  coordinates  of  the  source. 

Therefore 

W  (r,r',  W;  KR1  -  L  (ui>rL6lI  G  (KR>Z)  exp  (j™  -KRr"-  KRr  cos  (0-0") 

(r')l/2  K  4 

4.12 


Now  let  us  consider  a  continuous  isotropic  source  distribution  in  the 
interval  a  <  r'  <  b,  Wr  then  may  be  expressed  as 


/2  ^  ^ 

exp  ^-jKRrcos(0-6 '))  dO'y  (r")1/2L(w,r')  exp  (-jKRr')dr' 


=  G  (KR,Z)ejJJ0  (KRr)  L  (u,,KR;a,b) 


4.13 


where 


/b 

(r')1/2  L  (u>,r ')  exp  (-jKRr')  dr' 


4.14 


Now  let  us  assume  that  the  vertical  strain  field  is  observed  at  the  point 
defined  by  the  radius  vector  ?  +  6 r.  Then  by  analogy  with  the  previous 
arguments 


6Wr 

6Z 


6G  (Kb,Z) 
6Z 


2-  J, 


(Kr| r  +  5r| )  L  (w,KR;  a,b) 


4.15 
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If  4>ei  is  the  cross  power  spectral  density  relating  the  strain  observations 
at  f  +  6?  to  the  displacement  observations  at  r  then 


4>ei  =  WR*  4.16 

6z 

=  |l|2  =  G  *  *£  J0  (KRr)  Jo  (Kr  | r  +  firj )  4.17 

6z 

Now  moving  the  vertical  displacement  point  to  the  origin  yields 

*ei  =  |  L| 2  G*  iG-  Jq(K  fir)  4.18 

6z 


By  similar  arguments  the  power  spectral  density  of  the  displacements, 
and  the  power  spectral  density  of  the  strains  are  found  to  be 


rn 


=  L 


4.19 


<l> 

ee 


4.20 
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5. 


NOISE  REDUCTION  IN  MULTICOMPONENT  FIELDS 


5.1  SCATTERED  RAYLEIGH  WAVES  CONTAMINATED  BY  STRAIN  SYSTEM  NOISE 

Let  R  be  the  ratio  of  the  noise  power  at  the  output  of  the  prediction  error 
processor  relative  to  the  noise  at  the  output  of  the  inertial  sensor.  Then 
according  to  this  definition  complete  noise  rejection  is  achieved  when  R  =  0 
and  increasingly  greater  amounts  of  noise  are  passed  by  the  processor  as 
R  ■*  1.  From  equation  1.6 


Let  us  assume  now  that  the  earth  noise  consists  entirely  of  Rayleigh  waves 
from  scattered  surface  sources  and  that  the  data  observed  at  the  output  of 
the  strain  sensor  are  contaminated  by  system  noise.  If  the  strain  and 
intertial  sensors  are  colocated,  then  from  equations  4.17  through  4.19  relevant 
power  and  cross  power  spectral  densities  are  given  by 


li 


LM  G 


5.2 


4>ee  -  I  L 1  2  |G  I2  ♦  | Ng | 2 


rei 


|l|2  g*  «£ 

6Z 


5.3 

5.4 


where  |Ng|2  is  the  system  noise  power  at  the  output  of  the  strain  sensor. 
Therefore  in  this  case  is  given  by 


1 

1+A 


5.5 

5.6 
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Let  B  define  the  ratio  of  system  noise  to  total  noise  power  at  the  output 
of  the  strain  sensor.  Then 


and 


A  = 


B 

1-B 


Therefore 


Y2 . =  1-B 

ei 


5.7 


5.8 


5.9 


R  =  B 

In  other  words  in  this  particular  case  the  reduction  in  noise  achieved  through 
prediction  error  filtering  using  a  linear  combination  of  the  outputs  of  strain 
and  inertial  sensors  is  limited  only  by  the  presence  of  system  noise  at  the 
strain  output. 


5.2  SCATTERED  RAYLEIGH  WAVES  CONTAMINATED  BY  WIND  NOISE 


In  this  case  let  the  power  spectra  of  the  Rayleigh  wave  noise  at  the  output 
of  the  inertial  and  strain  sensors  be  denoted  by  |Wp|2  and  |Wn|2  respectively. 


Similarly  let  |Ww|2  and  |W^|2  be  the  power  spectra  of  the  wind  noise  at  the 


outputs  of  the  inertial  and  strain  sensors, 
ratios  Q  and  T  by 


Furthermore  let  us  define  the 


Q2  - 


K!2 

Iwwl2 


5.11 


T2  = 


1  wr1  2 
|"rI2 


From  arguments  presented  by  Sorrells  (1971)  Q2  can  be  shown  to  be  given  by 


Q2 


5.12 
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where  H  is  the  quasi  static  vertical  displacement  response  of  the  earth  to 
atmospheric  pressure  variations  and  is  the  vertical  strain  response.  For 

6  Z 

observations  at  the  surface  of  a  perfectly  elastic,  homogeneous  and  isotropic 
half  space 


H  =  A 

c0  f X+2yA 

5.13 

2p 

to  V.X+y  J 

5H  _  1 

(  M  ^ 

6Z  2  y 

\,X+y  J 

where  C  is  the  local  wind  speed.  Therefore 

2 

n 

Q2  -  4 

0T 

0® 

5.14 

From  equations  3.34 


♦  2Kr 

”'2Kr  -  k6>  -  2Kr 


S.  15 


and  for  the  case  where  X  =  y 


T2 


16  B2 

to2 


5.16 


Now  using  equations  5. 11  to  express  the  inertial  spectral  components  in  terms 
of  their  strain  counterparts  it  is  easy  to  show  that  in  this  case  Y2,is  given 
by 


(Qlw*l2  ♦  Tl wr!2) 

(q|W'|2  ♦  T|W'|2)2+  (Q-T)|W'|2|Wr|2 


5.17 
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Then  letting  B  be  the  ratio  of  wind  noise  power  to  total  noise  power  at  the 
output  of  the  strain  sensor;  i.e 


l«wl2 


B  = 


|w;|2  + 


it  can  be  shown  that 


R  =  1  -  y2 
ei 


(>-?) 

|2B(1-B) 

2  B(l- 

B)  * 

V  (l-§  BJ2 

and  from  equations  S.14  and  S.16 


a  =  _jl_ 

T  16p 

Now  c  is  generally  of  the  order  of  meters/sec  while  3  is  of  the  order  of 
kilometers/sec.  Thus 


=  0  CIO*3) 


Therefore 


1  -S 

T 


,  1 


5.22 


and  for  B  <  1 


5.23 
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A 


In  othes  words  for  small  values  of  B,  the  addition  of  wind  noise  to  the 
Rayleigh  wave  field  is  comparable  to  the  addition  strain  system  noise,  insofar 
as  its  impact  on  the  noise  reducing  capabilities  of  a  two  channel  strain 
inertial  processing  system  is  concerned.  Notice  that  this  approximation  is 
no  longer  valid  for  values  of  B  near  1,  i.e.,  when  the  strain  noise  is 
completely  saturated  by  the  wind  induced  component.  Under  these  conditions 
it  is  convenient  to  write  R  as 


R 


r/  (1  -  SB)  SB 
C  (1-SB)  SB  +  1 1  -  r.  ( 1  -  SB)  ] 2 


who  re 

,  .  i  -  a 

r 

and 

B  =  1  -  SB 


Assuming  SB  <<  1  then  equation  5.J4  can  be  approximated  by 


R  = 


which  shows  that  R 
Summarizing  then 


0  as  SB  *  0  for  2-/0 

T 

R  '  B  B  -  1 
R  ►  0  as  B  *■  0 


and 


R  : 


R  *  0  as  B  >  I 
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it 


Equations  5.26  and  5.27  state  that  the  noise  rejection  provided  by  a  strain 
inertial  system  will  be  complete  only  if  the  noise  at  the  output  of  the  strain 
transducer  consists  of  either  pure  Rayleigh  wave  noise  or  pure  wind  generated 
noise.  Otherwise  the  noise  reducing  capabilities  of  the  system  will  be 
functinally  dependent  on  the  percentage  of  wind  noise  in  the  noise  observed 
at  the  output  of  the  strain  sensor. 

Finally  it  can  be  readily  shown  that  R  has  a  maximum  value  at 


at  this  point 


5.28 


5.29 


For  the  case  of  Rayleigh  waves  contaminated  by  wind  noise  the  maximum  value 
of  R  occurs  at  B  ;  1  and  at  this  point  R  max  ~  1.  Thus  except  for  a  very 
small  range  of  values  near  B  =  1 

R  =  B  0  v  B  <  1  5.50 

which  is  identical  to  the  relationship  found  for  the  strain  system  noise. 


5.5  SCATTERED  RAYLEIGH  WAVES  CONTAMINATED  BY  P  WAVES 

An  expression  similar  to  5.19  may  be  derived  for  the  case  where  the  propagating 
earth  noise  consists  of  both  scattered  P  and  Rayleigh  waves.  In  this  instance, 
however,  it  is  instructive  to  express  the  noise  reduction  in  terms  of  the  ratio 
of  the  P  wave  noise  power  to  the  total  noise  power  observed  at  the  output  of 
the  inertial  sensor. 

Let  Bp  denote  this  ratio.  Then 

B  -  - !!tll - -  5.3] 

*  I«rP 

where  |Wp|?  is  the  P  wave  spectral  noise  power. 
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Let  the  parameters  P  and  V  be  defined  by 


|p|2  ,  P2 


V|2  u  \2  . 

|wRl2 


5.32 


5.33 


where,  as  before,  the  primed  quantities  denote  spectral  densities  at  the 
output  of  the  strain  sensor. 

Then  it  can  be  easily  shown  that 


R 


\ '  *0]  »r  (1-Bp) 

1  -[•-  (ff]"r 


5.34 


Now  suppose  we  confine  our  attention  to  the  case  where  the  observations  arc 
made  at  the  surface  and  the  scattered  P  waves  arrive  at  the  site  with  near 
vertical  angles  of  incidence.  This  approximates  the  so  called  "mantle  P  wave" 
noise  field.  Then 


P|- 


r2  u  ♦  n  ♦  Q 


V1  -  n  -  a 


5.35 


where  F  ;md  si  are  defined  by  equations  3.10  through  3.13 
and  assuming  X  =  u. 


ip)  2  ;  4  u>‘  s  i  n  (  3  sin- 0-1  I 

9  a2  }  (3  sin2e-)2  f 


5 . 36 
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where  0  is  the  angle  of  incidence  for  the  P  wave  noise.  Then  using  equation 
5.16 


I  2  =  — 


1  .  0>2 
T2  '  1662 


5.37 


and 


P|2  P2  16  sin2e  (3  sin2  e-1 


I  Vi2  V2 


27 


sin2  e-l): 


5.38 


If  we  confine  our  attention  to  near  vertical  angles  of  incidence  n  <<|and 
equation  5.37  reduces  to 


fl  =  16  e2 

V2  27 


and  given  this  condition  R  reduces  to 


R  =  Bp  ,  Bp  t  1 


Near  Bp  =  1  we  may  replace  Bp  by  l-6Bp  and 
neglecting  power  of  <5Bp  >  2  one  obtains 


5.39 


5.40 


5.41 


It  should  be  noted  that  according 


to  the  definitions  of  Bp  and  SBp 


=  _ 1^1  2  ... 

|Wp|2  ♦  | WR | 2 

which  is  the  fractional  Rayleigh  wave  power  at  the  output  of  the  inertial 
sensor.  Let  this  ratio  be  defined  by  BR.  Then  equation  5.41  may  be  rewritten 

r  :  |l+  (£)  W  Br  0  5.42 
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It  can  be  easily  shown  for  this  case  that  R  has  a  maximum  at 


BP 


1 


1 


+ 


P 

V 


and  at  this  point 


Since  —  <<  1 ,  R  reduces  to 
V 


R  *  1  -  ifi 

|V| 


5.44 


5.45 


Thus  where  the  propagating  major  contains  both  Rayleigh  wave  and  mantle  P  wave 
components,  the  strain  inertial  processing  system  will  completely  suppress 

I  P  i 

the  Rayleigh  wave  component  in  the  interval  0<  Bp<  1  -  J _ L.  As  B.>  increases 

| p |  lvl 

through  1  -  j-y-j  the  system  will  begin  to  suppress  the  P  wave  component  and 
total  noise  rejection  is  achieved  at  Bp  1. 
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